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A theory of critical fluctuations in extreme type-II superconductors subjected to a finite but weak 
external magnetic field is presented. It is shown that the standard Ginzburg-Landau representation 
of this problem can be recast, with help of a novel mapping, as a theory of a new "superconductor", 
in an effective magnetic field whose overall value is zero, consisting of the original uniform field and a 
set of neutralizing unit fluxes attached to iV$ fluctuating vortex lines. The long distance behavior of 
this theory is governed by a phase transition line in the (H,T) plane, T<s>(H), along which the new 
"superconducting" order parameter, $(r), attains long range order. Physically, this novel phase 
transition arises through the proliferation, or "expansion", of thermally-generated infinite vortex 
loops in the background of field-induced vortex lines. Simultaneously, the field-induced vortex lines 
loose their effective line tension relative to the field direction. It is suggested that the critical behavior 
at T<s,(H) belongs to the universality class of the anisotropic Higgs-Abelian gauge theory, with the 
original magnetic field playing the role of "charge" in this fictitious "electrodynamics". At zero field, 
$(r) and the familiar superconducting order parameter, ^(r), are equivalent, and the effective line 
tension of large loops and the helicity modulus vanish simultaneously, at T = T c o. In a finite field, 
however, these two forms of "superconducting" order are not the same and the "superconducting" 
transition is generally split in two branches: the helicity modulus typically vanishes at the vortex 
lattice melting line, T m (H), while the line tension and associated "I>-order disappear only at T&(H). 
We expect T$(H) > T m (H) at lower fields and 2* (if) = T m (H) for higher fields. Both $- and 
^-order are present in the Abrikosov vortex lattice (T < T m (H)) while both are absent in the true 
normal state (T > 7$ (H)). The intermediate ^-ordered phase, between T m (H) and T§{H), contains 
precisely N§ field-induced vortices having a finite line tension relative to H and could be viewed as 
a "line liquid" in the long wavelength limit. The consequences of this "gauge theory" scenario for 
the critical behavior in high temperature and other extreme type-II superconductors are explored 
in detail, with particular emphasis on the questions of 3D XY versus Landau level scaling, physical 
nature of the vortex "line liquid" and the true normal state, and fluctuation thermodynamics and 
transport. It is suggested that the empirically established "decoupling transition" may be associated 
with the loss of integrity of field-induced vortex lines as their effective line tension disappears at 
T<s>(H). A "minimal" set of requirements for the theory of vortex lattice melting in the critical 
region is also proposed and discussed. The mean-field based description of the melting transition, 
containing only field-induced London vortices, is shown to be in violation of such requirements. 

Pacs-numbers: 74.40. +k, 74.25.Bt, 74.20.De, 74.25.Dw, 74.25.Ha, 74.60.Ec 



I. INTRODUCTION 

Recent intense activity in the area of superconducting 
fluctuations has brought into sharp focus the following 
fundamental questions: .-What is the relationship between 
the Landau level-basedW and the 3D XY-basedHTJ de- 
scriptions of superconducting fluctuations in magnetic 
field? Can the mean-field based London model contain- 
ing only magnetic field-induced vorticesQ describe the 
vortex lattice melting transition in the region of strong 
(critical) fluctuations? What is the nature of the nor- 
mal phase and can it be usefully represented as a "line 
liquid" cl of field- induced vortices? What role is pl£Wf4 
at finite fields by thermally-generated vortex loopsJ§t2l 
which are responsible for the zero-field transition in ex- 
treme type-II superconductors? Particular importance 
and urgency has been attached to these questions follow- 



ing the ground-breaking experimentsUH3 on the ther- 
modynamics of vortex lattice melting transitions which 
clearly indicate that the low-field end of the melting line 
is entering the critical regime of high temperature super- 
conductors. 

In this paper, precise answers to these questions are 
provided within a theoretical framework which allows for 
a systematic solution to the problem of critical fluctua- 
tions in an extreme type-II superconductor subjected to a 
finite, but weak magnetic field. This framework is built, 
around the "gauge theory" scenario proposed earlier.Ll 
Two main predictions follow from this scenarios first, 
there is a new transition line in the H — T phase di- 
agram, T$(H), along which a thermally-generated vor- 
tex loop "expansion" takes place, reminiscent of the zero 
field transition. At T$(iJ), well defined field-induced vor- 
tex lines are formed, having a finite line tension relative 
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to the field direction. Initially, these lines are in a liq- 
uid state and solidify only at some lower temperature 
T m (H) (Fig. 1). This is different from the Abrikosov's 
theory, where vortices and their (Abrikosov) lattice are 
formed simultaneously, at H C 2(T); second, in contrast to 
the 3D XY behavior at zero field, the description of the 
critical behavior along T$(H) requires the combination 
of a complex "superconducting" order parameter $ as- 
sociated with vortex loops and a fictitious gauge field S, 
describing fluctuations in the background system of field- 
induced vortex lines. The magnetic field determines the 
"charge" which couples $ and S. The physical picture 
arising from the "gauge theory" is remarkably detailed 
and compelling—and so entirely distinct from the "stan- 
dard" approachEJ, that a concentrated effort should be 
directed at exploring its consequences. The main pur- 
pose of this paper is to provide a novel explicit model for 
critical fluctuations, to examine its main ramifications 
in some detail, and to advance a set of specific predic- 
tions which can help establish the value of the "gauge 
theory" description though experiments and numerical 
simulations. 

The essential feature of our description is that it con- 
tains, on equal footing, both the field-induced vortex lines 
and thermally-generated critical fluctuations of the su- 
perconducting order parameter ^, associated primarily 
with vortex loops. It is theJatter that dominate the en- 
tropy in the critical region.EI This is in fundamental con- 
trast to other approaches which include only the S^lfLi 
induced vortices: the Landau level (LL) descriptionE£Tt£l 
(where other fluctuations become irrelevant at suffi- 
ciently high fields,) and the mean-field based picture of 
London vorticest3 (where other fluctuations can be ig- 
nored at sufficiently low temperatures, far below critical 
regime). Following the prescription proposed earliercl, 
which seeks to conveniently isolate the background of 
field-induced from thermally-generated degrees of free- 
dom, we derive the following results: First, it is shown 
in Sec. II that the familiar and frequently used "helium" 
or "London model" of extreme type-II superconductors, 
in which the amplitude fluctuations are suppressed, al- 
lows for a direct mapping of the original problem to that 
of a new "superconductor", whose order parameter $ 
experiences an overall magnetic field composed of the 
uniform external field, H, and the set of iV$ neutraliz- 
ing "fluxes" attached to fluctuating vortex lines. This 
novel mapping constitutes an explicit and transparent 
realization of the general connection proposed in Ref.tl 
The "helium model" is then a candidate to, in addition 
to the familiar vortex lattice melting line, T m (H), ex- 
hibit the conjectured "^-transition" (or the vortex loop 
"expansion" transition in a finite field), the universal- 
ity class of which is djefined by an anisotropic Higgs- 
Abelian gauge theoryflcj Physically, this ^-transition 
corresponds to vanishing of the effective line tension for 
very large thermally-generated vortex loops: at T$(iJ), 
the energy-entropy balance in the free energy shifts in 
favor of large loops and spontaneously created infinite 



vortex-antivortex paths proliferate across the system. 
The ensuing change in the topology of the vortex paths 
results in the breaking of a "dual" U(l) symmetry and 
a thermodynamic phase transition. Simultaneously, the 
field-induced vortex lines loose their line tension relative 
to the field direction and the "line liquid" description 
breaks down. As our second result, it is shown that the 
fictitious gauge theory passes a crucial test, allowing us 
to connect its "charge" to the original magnetic field. 
Third, we use this connection in Sec. IV to construct 
scaling functions for the critical thermodynamics of ex- 
treme type-II superconductors. Furthermore, the much- 
debated difference between the 3D XY-like description 
at low fields and the LL description appropriate at high 
fields is closely linked here to the difference between the 
extreme "type-II" and the extreme "typeT" behavior of 
the gauge theory (Sec. III). Fourth, it is shown in Sec. V 
that a vortex loop "expansion" leads to an abrupt drop 
in the coefficient of the q 2 -term in the helicity modulus, 
from which one can extract the thermodynamic exponent 
{v) of the ^-transition. Related criteria are also proposed 
which test for the presence/absence of effective "diffu- 
sion" of vortex lines along the fidd and demonstrate close 
relation betweenr-the "$-order"B and viability of the vor- 
tex "line liquid" □ description. These predictions, based 
only on global topological properties of loops and lines, 
can be used to efficiently identify the vortex loop "expan- 
sion" line, T$ ( H ) , in numerical simulations of the weakly 
frustrated 3D XY and related models. Fifth, assuming 
widely used form of dynamical scaling, the explicit ex- 
pression for the fluctuation conductivity, a(T,H), is de- 
rived in Sec. VI in the vicinity of the T${H) line. At 
T$(H) there is an experimentally detectable rapid onset 
of additional dissipation, caused by thermally expanding 
vortex loops whose size is reaching sample boundaries. It 
is tempting to associate this onset at T$ (H ) with what 
is empirically known as the "decoupling" transition,!^ al- 
though the physical origin of such additional dissipation 
in our theory is-entirely unrelated to any "decoupling" 
of any "layers" £3 Instead, it signifies the loss of integrity 
of field-induced vortex lines as their effective line tension 
disappears at T$(iJ). At this point, one also expects 
a distinct change in the pinning properties of the liq- 
uid state: there is no pinning in the true normal state 
above T$(H) just like there is no pinning above T c q. 
In this sense T$(H) represents an upper boundary for 
pinning and could be viewed loosely as "renormalized" 
H C 2{T). Finally, in Sec. VII, it is demonstrated that 
the vortex lattice melting transition in the critical region 
involves simultaneous ordering of the field-induced and 
thermally-generated degrees of freedom and thus cannot 
be faithfully represented by a mean-field based London 
modeljQ which includes only the former. Actually, as the 
melting line tends toward T c q in the limit of vanishing 
magnetic field, the entropy change involved in ordering 
of thermally-generated loops overwhelms the configura- 
tional entropy of the field-induced vortex lines. This 
provides direct theoretical support for the fundamental 
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nature and significance of the experiments, by Zeldov et 
al.Jld Schilling et al.J13 and Roulin et al., Ill and new nu- 
merical simulations of Nguyen and Sudbd^J'EJ. 

II. FROM THE GINZBURG-LANDAU THEORY 
TO THE GAUGE THEORY 

The starting point is the anisotropic Ginzburg-Landau 
(GL) theory Z = / £>* exp{- / d 3 rF/T}, where 

T = a\^ + £ 7mI(V m + ^A m )*| 2 + ^| 4 , (1) 

/Hl.-L 

and a = oq(T — T c ), 7^, and (3 are the GL coefficients. 
Free (periodic) boundary conditions are imposed in || (_L) 
direction. The limit k — > 00 is considered, which is par- 
ticularly appropriate for high temperature superconduc- 
tors (HTS). In this limit, the external magnetic field, 
H = V x , acts as a constraint, forcing every allowed 
configuration of the system to have the overall vorticity 
7V$ along H (|| z). The overall vorticity along H is de- 
fined as a line integral J dl ■ V<p/2ir, where the contour 
of integration goes around perimeter of the system in the 
xy-plane and ip(r) is the phase of N^>, the number of 
elementary flux quanta <fio, is given by L? L /2ir£ 2 , where 
I = ^c/2\e\H is the magnetic length. It is assumed that 
this constraint is enforced by N$ vortex paths, meander- 
ing from one end of the system to another, along H.lI 

In this paper, a new method is introduced 
to enforce the constraint explicitly, by con- 
sidering a different partition function: Z' — 
fV$J nil* (Pn [a]/JV«!) exp{- / d^rP/T) , with 

lei B 
T' = a|$| 2 + 7ai |(V m + zU, + — A M )$| 2 + ||$| 4 . (2) 

Z' describes the system of iV$ "shadow" , or s, vortices 
{ri[s]} in thermal equilibrium with a complex field $(r). 
These s vortices sample arbitrary paths that originate 
(terminate) at z = (z — L\\) and differ from the ones 
introduced in RefE by the full inclusion of "overhang" 
configurations. The effective magnetic field H' expe- 
rienced by $ consists of the uniform external field H 
and the collection of unit "fluxes" attached to s vortices: 
V x U = 27rn s (r); V-U = 0, where n s (r) is the flux den- 
sity associated with a given configuration of s vortices, 

{ri[*]}£a 

n a (r) = / dvi5{r - Yi[si]) , (3) 
i J £ 

with C denoting the line integral. The net value of H' 
averaged over the system vanishes. 

The superconductors (Q) and (|2|) are equivalent within 
the familiar "helium model" of extreme type-II behavior; 
they are just two different representations of the same 



physical problem. To show this we recall the main fea- 
tures of the "helium model" J?J the true transition tem- 
perature T c q (Fig. 1) is assumed to be sufficiently below 
the mean-field T c for amplitude fluctuations to have effec- 
tively subsided. Around T c q, the relevant fluctuations are 
considered to be those of London- type vortex loops/lines 
with steric repulsion and well-defined, tight cores of size 
a <ti £. Consider now a single configuration of these loops 
and lines. First, we extract the singular part of Vy(r) by 
solving two equations: Vx Vi/J = 27rn(r) and V- V92 = 0, 
where n(r) is defined by the same expression as n s (^) but 
with the summation running over all vortex loops and 
lines. After this "vortex" part has been extracted, the 
rest of 'J'(r) is assumed to take the form which minimizes 
J- for a given configuration of these line singularities. We 
then integrate over all regular ( "spin- wave" ) fluctuations 
in <p(r). Finally, all such distinct configurations of vor- 
tex loops and lines are summed over to produce Z (Q). 
This same procedure is imposed on $(r): first we extract 
the part of its phase, V0(r), due to vortex loop/line sin- 
gularities in $(r) and then determine the rest of 3>(r) 
by minimizing T' (|2|) for a given configuration of these 
line defects and s vortices. Again, we integrate over all 
"spin-wave" fluctuations of <fi relative to this given config- 
uration of loops and lines. By direct comparison of these 
"helium model" expressions obtained from and (j^), 
it is evident that all configurations contributing to the 
original Z are reproduced in Z' and have the same en- 
ergy. However, some of these configurations are counted 
more than once in Z' . This overcounting of configura- 
tions in Z' relative to Z, given by (A$ + N a )\/N^\N a \ 
with 7V a being the number of vortex lines in $ which tra- 
verse the sample along H, is a surface effect in 3D and 
should be unimportant in the thermodynamic limit, L±, 
Lu — > 00. Moreover, within the conjectured <I>-ordered 
phase (Fig. 1), the configurations with N a ^ are irrel- 
evant in the thermodynamic limit and there is no over- 
counting at all. We conclude that, within the "helium" 
model, the free energy evaluated from Z' coincides with 
the free energy of the original problem (Q) and the two 
superconductors have identical thermodynamicsEj Con- 
sequently, Eq. (||) accomplishes a straightforward and 
transparent reformulation of the original problem, in the 
spirit of Ref.t], while avoiding more cumbersome gauge 
transformation method employed there.E£l More details 
on the "helium model" are presented in Appendix A. 

If we relax the above minimization condition on the 
amplitude of our order parameters, \& and <!>, and permit 
weak amplitude fluctuations, we expect that the above 
close relation between Z (Q) and Z' (||) still holds, as 
long as the parameters of the GL theory keep us in the 
extreme type-II limit. This requires the average core size 
a to be smaller than the average spacing between vortex 
segments, so that vortex excitations remain well-defined. 
It is precisely this same requirement that is invoked to 
justify the frequent use of the "helium model" to emulate 
fluctuation behavior of extreme type-II superconductors 
in zero field. It is natural to expect that, if such re- 
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quirement is satisfied at zero field, it will remain so at 
low fields, such that a <C £■ Based on this, on the equiva- 
lence of representations (|l|) and (||) in the "helium model" 
limit, and on our general expectation that the extreme 
type-II behavior with only weak amplitude fluctuations is 
effectively equivalent to the "helium model" , for the rest 
of this paper I consider (||) to be simply an alternative 
formulation of the original problem. This new reformu- 
lation (J2j) can now be used instead of (|l]) to compute 
various fluctuation properties and, most importantly, its 
critical behavior should coincide with that of the original 
GL theory (0). 

The advantage of Z' ^ is that, by isolating the back- 
ground of field-induced degrees of freedom (s vortices), 
it focuses our attention on the new "superconducting" 
order parameter, $(r), and its spatial correlations, mea- 
sured by ($(r)$*(r')), where (• • •) denotes thermal aver- 
age over Z' . All excitations of $(r) are thermally gener- 
ated, in the following precise sense: every configuration 
o/<I>(r), contributing a finite weight to Z' in the thermo- 
dynamic limit, has the overall vorticity along H equal to 
zero. In particular, $(r) contains vortex loop excitations, 
whose "expansion" across the system is the mechanism 
behind the H = superconducting transition (Fig. 2). 
By focusing on we can fashion a theory of the strongly 
interacting Wilson-Fisher. (3D XY) critical point, "per- 
turbed" by a weak ficldjj This is precisely the opposite 
of the classic approach,c3 where the Gaussian theory in 
a finite field is perturbed by weak interaction. Such ap- 
proach starts with the LL structure from the outset and 
its critical behavior is always dominated by the lowest 
LLEj'Ej In the new formulation (g) we had built in from 
the start our expectation that the weak field modifies 
zero field configurations only by introducing a low den- 
sity of s vortex lines, the cores of which are well-defined 
by virtue of strong amplitude correlations at the 3D XY 
critical point. This is a "low- field" approach by design 
and offers a better prospect of constructing the desired 
theory. 

To extract such a theory from (g) we must resort to ap- 
proximations. We construct the long wavelength (^> £) 
limit of (0) by coarse-graining vorticity fluctuations pro- 
duced by s vortices. The "hydrodynamic" vorticity V(r) 
is defined as the coarse-grained version of the "micro- 
scopic" flux density An s (r) = n s (r) — (2tt£ 2 )~ 1 z. Upon 
inserting J DV4[V(r) - An s (r)] in (g), integrating over 
{r,[s]}, and after introducing the fictitious vector poten- 
tial V x S = 27rV; V ■ S = 0, the effective long wavelength 
theory becomes ~ / X>$ / X>Sexp{- f d 3 rF cS /T}, 

^ eff = a |$|2 +7/i | jD/i $|2 + || $ |4 + ^ i(VxS) 2 ^ (4) 

where = V M + iS^ and K±(T,H) = c±T~ 1 T£, 
Ku(T,H) — c»TT£. Higher powers and derivatives of 
(V x S)„, essential for the description of the vortex lat- 
tice melting, also appear in (^), but are unimportant at 
T$(H). Detailed derivation of the "gauge theory" (Q) is 



given in Appendix A. Cjii ~ 0{\) are dimensionless and 
have a relatively weak H, T dependence in that portion, 
of the critical region which is well described by Eq. (||) 
A close relation between cmi and the components 

of the helicity modulus tensor of the GL theory (0) is 
discussed later in the text (see Sec. V and the Appendix 
B). T is the anisotropy at T c q. Small independent cor- 
rections to GL coefficients that also should appear in ([|) 
are ignored, since they are not important for our present 
purposes. 

The following assumptions have been used in going 
from (||) to (||) (see also Appendix A): 

i) The correlation length, £$, associated with the new 
order parameter <3>, is not limited by t and can be 
much longer than the original superconducting correla- 
tion length, £ sc , associated with *S>. Of course, this is the 
basic reason why we are interested in the reformulation 
(g) in the first place. When ^rfi £sc, this assump- 
tion enables us to drop as irrelevant!] at long distances, 
terms containing higher derivatives and powers of S from 
Eq. (|J). Note, however, that such higher order terms in 
(Q) , particularly those reflecting the absence of up-down 
symmetry along H ((V x S)| and the like), must be re- 
stored when discussing vortex lattice melting (Sec. VII). 
The gauge theory (|]) offers in this case (£$ 3> £, £ sc ) a 
direct access to the deeply non-perturbative regime of the 
original GL theory (|l|), characterized by £/£ ^> 1, where 
£ is the H = correlation length. In the opposite case, 

<C £, we are in the perturbative regime, £/£ -C 1, of 
the original theory. The long wavelength expansion that 
led from (|^) to (Q) is then not justified and the new re- 
formulation (g) is not particularly useful. 

ii) The system (|l|) is not in its superconducting state 
in the vicinity of the putative ^-transition. This as- 
sumption fixes the form of the last two terms in (Q) 
(see Sec. V and the Appendix B). Physically, it means 
that the system of s vortices contains configurations that 
wind from one end of a sample to another in the _L di- 
rections (xy-plane). The presence of such windings al- 
lows complete "screening" of an arbitrary infinitesimal 
field h(r) added to H and the helicity modulus ten- 
sor vanishes along all of its principal axes (see Sec. V 
and the Appendix B). | _This assumption has a strong 
theoretical justification.^] It must be emphasized, how- 
ever, that, to my knowledge, there is no rigorous ar- 
gument which could rule out another possibility, that 
of the state right below the ^-transition being an ex- 
tremely anisotropic "superconducting" liquid, contain- 
ing no windings in the xy plane, with a finite helic- 
ity modulus along the field and zero perpendicular to 
it. Indeed, some numerical studies are suggestive of 
this second—^ssibility.EJa However, other numerical 
simulationsoOO, as well as the available experimen- 
tal data, favor our original non-superconducting liquid 
assumption. Both alternatives can be described within 
the framework of the gauge theory, with K±/Ku — > oo 
and potentially finite "mass terms" below T$(H) added 
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to (Q) in tha-epctreme anisotropy case. On general physi- 
cal grounds,E3 I have chosen to explore in this paper the 
case of finite anisotropy ratio K±/Kn and vanishing mass 
terms but the reader should be aware that-the extreme 
anisotropy alternative remains a possibility^ and would 
lead to results which, while similar on general level, differ 
in details.E-3 

iii) There are two relevant lengthscales controlling the 
critical behavior: £$, which characterizes both the spa- 
tial correlations of $ and the size of "overhangs" in the 
system of s vortices, and £, which characterizes the long 
wavelength fluctuations of the background field-induced 
vorticity. 

iv) The core effects can be ignored. Clearly, the "he- 
lium model" itself is perfectly well defined in the limit 
a — > 0. For a small but finite, there is a small correc- 
tion to the core line energy e c — » e c + w c H • v, where 
v = dr/ds is the "velocity" of a vortex segment and 
w c H/e c ~ a 2 /£ 2 -C 1. Similarly, there are "velocity" de- 
pendent corrections to the short range repulsion between 
vortex cores. Such terms are irrelevant since they result 
in higher order derivatives in (0). For example, it is easy 
to see that the correction to the core energy cancels out 
for any finite vortex loop and can be factored out for s 
vortices. 

What is the physics behind gauge theory (||)? The 
external field has been eliminated from the gradient 
terms in @) ((H') = 0) but, of course, it has not 
vanished: it reappears through the TJ-dependence of 
K±n. The gauge theory (Q) can be viewed as ficti- 
tious, anisotropic "electrodynamics" with "magnetic per- 
meability" no = \/AttT. The "vector potential" S is cou- 
pled to the "matter" field $ via "electrical charge" 



pl/3 



cx Vh 



The above "charge" and K±,\\ describe the "polarizabil- 
ity" of the medium composed of s vortices and are di- 
rectly related to the long wavelength components of the 
helicity modulus tensor (Sec. V and the Appendix B). 
This picture embodies the physical idea that the dom- 
inant effect of weak magnetic field in ([!]), once £ sc has 
saturated to ~ £, arises through "screening" of large 
thermally-generated loops by the background of field- 
induced vorticity, at distances 3> £■ Such "screening" 
reduces the effective line tension of these large loops rel- 
ative to its value at the H = (e± t \\ = 0) transition. The 
strength of the "screening" is measured by the fictitious 
"Ginzburg parameter" of (0), 



/3£ 



1 



H 



(5) 



where c = {c 2 ^) 1 ^ and £ GL = (£gl±£gl\\) 1/3 , with 



£gl_l,|| = \/l±,\\/0'oTc being the GL coherence lengths. 
b — P/a 2 ) T c 3L and q 2 = e 2 £_GL are the dimensionless 
quartic coupling and "charge" , respectively. As H — > 0, 



the fictitious "charge" vanishes and we recover the zero- 
field 3D XY critical point. For H finite but weak, the 
"screening" is weak (k s 3> 1), indicating that the effects 
of finite e are small compared to strong amplitude cor- 
relations produced by the quartic term in the GL theory 
(1). We have therefore manufactured a critical theory 
( 1 ) describing the strongly- interacting Wilson-Fisher (3D 
XY) critical point weakly "pertpjbed" by a finite mag- 
netic field (finite "charge" e_u|)£3 

The gauge theory scenario is clearly different from 
what takes place in spin systems, where the external 
field couples paramagnetically to the order parameter. In 
an extreme type-II superconductor (jl|), the diamagnetic 
coupling of H to f does not explicitly break the U(l) 
symmetry, which was spontaneously broken at the zero- 
field 3D XY critical point. The high temperature phase 
(true normal state) still retains the full U(l) symmetry. 
This symmetry can be broken at low temperatures, ei- 
ther in a "simple" way, with \l/ acting as the order pa- 
rameter, as is the case in the "vortex solid" state, or 
in a more subtle fashion, with $ assuming the role of 
the new order parameter. Similarly, the gauge theory 
(|j) differs from frequently used "dimensional reduction" 
approaches ,E3 where the behavior of ([!]) at finite fields is 
related to that at zero field but in a finite system, the 
size of which is set by the magnetic length, £. A typi- 
cal dimensional reduction (D — > D — 2) approach leads 
to the superconducting correlation length which is lim- 
ited by £, i.e. the "svstem size" . This agrees with the 
gauge theory scenariojj since "electrodynamics" (Q) also 
predicts £ SC (H) ~ I in the critical region (see Sec. V). 
However, a dimensional reduction approach also predicts 
that all other correlations are limited by £, and, conse- 
quently, eliminates the possibility of any true thermo- 
dynamic phase transition in the GL theory (|lj). This is 
in contradiction with the overwhelming experimental and 
numerical evidence indicating some form of a "vortex liq- 
uid" to "vortex solid" transition at low temperatures. In 
sharp contrast, gauge theory (^) and reformulation (||) 
are fully three-dimensional theories, just like ([!]). They 
naturally lead to two basic types of correlations that can 
extend over distances 3> £ and produce phase transi- 
tion^) at low temperatures (Fig. 1): those associated 
with positional order of s vortices and the familiar super- 
conducting order parameter ^(r) and those associated 
with the new "sueerconducting" order parameter <I>(r). 

The conjecturea that connects the critical behavior of 
an extreme type-II superconductor (^) to a fictitious su- 
perconductor in zero field (0), the "charge" of which is 
set by the original external field H , must pass the fol- 
lowing test: the way H enters in must be consistent 
with its being a relevant operator of scaling dimension 2 
in the renormalization group (RG) sense at the 3D XY 
critical point. This-scaling dimension is suggested by di- 
mensional analysisEJ (relevant effects of H enter through 
the dimensionless ratio Q C LH = 0)/£ 2 cx H^ 2 C {H = 0)), 
is correct to two-loop ordertS and is likely an exact prop- 
erty of the original GL theory (|l]) by virtue of gauge 
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invariance. In addition, the scaling dimension appears 
independent on the nature of the zero-field critical point 
(i.e., whether it is 3D XY or Gaussian). On the other 
hand, as the "charge" e is turned on in the gauge the- 
ory (Q), the RG analysis indicates that first, the fpj|tG 
charge anisotropy (e_L ^ eii) is marginally irrelevant^B'c3 
and second, the scaling dimension of "charge" at the 3D 
XY critical point is 1/2, i.e., the relevant dimensionless 
operator is e^J ^ sc (e = 0). The second statement is exact 
to all orders in perturbative RG and is also independent 
on the nature of the neutral critical point.EJ Since, in 
T c q @), e 2 oc l/l oc \[H this translates immediately to 
the scaling dimension of H being 2, as required. More 
generally, for dimension D < 4, e 2 oc £ D_4 , while the 
scaling dimension of e is 2 — (Z)/2), again consistent with 
the scaling dimension of H being 2. Note that in both 
formulations, Eqs. (Q) and the corresponding rele- 
vant operators H and e (oc H 1 / 4 ) are protected against 
acquiring anomalous dimensions by the same symmetry, 
the gauge invariance. These results demonstrate internal 
consistency of the coarse-graining procedure leading to 
(||) and strongly support the conjectures that J- e g cap- 
tures the long wavelength (critical) behavior of (§) and 
(|l|). In what follows, I promote this conjecture to a fact 
and examine its consequences. 



III. HIGH FIELDS VERSUS LOW FIELDS 



Two key consequences for the physics of the present 
problem follow from Ted- First, the gauge theory ([|) has 
two distinct regimes of behavior: the weak "screening" 
limit (k s » 1) corresponding to the extreme "type-IP 
limit of the fictitious "electrodynamics" and the strong 
"screening" limit (k s C 1) corresponding to the extreme 
"type-I" behavior. The extreme "type-II" behavior of 
(@) isprecisely the low-field regime of the original the- 
ory (p which exhibits the 3D XY-like critical fluctua- 
tions. In this low-field regime, the "screening" provided 
by the background of field-induced vorticity is weak and 
the dominant fluctuations are still London-type vortex 
loops and lines. The core size a remains small and well- 
defined, kept in check by strong amplitude correlations 
coming from the quartic term in (|l|), just as was the 
case at the zero-field 3D XY critical point. It is in this 
sense (k s 3> 1) that we can think of a 3D-XY critical 
point weakly "perturbed" by a finite field.EJ In the ex- 
treme "type-I" limit, the situation is entirely different. 
There, the "screening" is strong (k s <C 1) and the am- 
plitude fluctuations ran rampant. It is not possible any 
longer to think of relevant fluctuations in the gauge the- 
ory ([I|), nor in (||) and (|l|), as being London-like vortices. 
Rather, amplitude fluctuations are now of essential im- 
portance and individual vortex cores are ill-defined. In 
the gauge theory (Q), the two regimes are separated by 
the condition k s ~ 1. However, a word of caution must 
be inserted here since, once we are in the "type-I" regime 



of , our original line of reasoning that lead from Eq. 
(|l|) to the gauge theory (|J) via reformulation @ , is itself 
compromised and it is not clear whether there is a useful 
connection between the extreme "type-I" limit of and 
our original problem ([j]). Instead, we must return back 
to the beginning ([j]) and start from scratch. It is natural 
to identify this extreme "type-I" behavior at high fields, 
characterized by strong amplitude fluctuations, as the 
regime in which the Landau level structure of the orig- 
inal GL theory (|l|) becomes important. The condition 
k s ~ 0.4/ \/2,E!J separating "type-II" from "type-I" "elec- 
trodynamics" in Eq. (Q), translates to the criterion for 
the external magnetic field, H ~ H s , telling us whether 
H is "low" or "high" . From Eq. (§) one g ets: 



^0.16 ; 



(0) 



(6) 



If H <C H s then the field is "low" and the use of a 3D 
XY-like description is justified. In the opposite limit, 
H > H s , the field is "high" and a 3D XY-like descrip- 
tion falls apart (Fig. 1). If this is the case, we must 
abandon our zero- and low-field imagery of the "helium 
model" and use as a starting point an approach that is 
explicitly designed to deal with a hi gh i fi eld behavior, an 
example being the GL-LLL theory.ll3i!3 Note that H s , 
within factors of order unity, coincides with Hi,, the field 
below which the high-field, Landau level-based descrip- 
tion breaks down, due to strong LL mixing.ES. Since this 
criterionll3 is derived from entirely different arguments, 
we briefly reproduce it here for completeness. Going back 
to the GL theory (|l|), we expand \I/(r) = J2j=o ^j( T ) m 
the set of LL manifolds, (r). Recast in terms of dimen- 
sionless variables, the GL free energy functional becomes: 



3=0 



(7) 



where t = (T/T c ) - 1, h = H/H% L {Q) and b is defined 
below Eq. (|5|). After rescaling \I/ and r by b in such a 
way that the coefficients of quartic and gradient terms in 
(0) are set to 1/2 and 1, respectively, the "mass term" 
for \E" j becomes 

t x h 

As we reduce the field H, for T in the critical region 
(T « T c ), the mixing of LLs becomes strong when h/b 2 
becomes some number of order unity. We can view this 
as a "Ginzburg criterion" along the iJ-axis. It implies 
that the high-field, Landau level description becomes in- 
adequate for fields less than: 



b 2 Hg L 



(0) 



(0) 



(8) 



where we have used a close relation between b apd a con- 
ventional Ginzburg fluctuation parameter, GiJla While 



G 



the definition and meaning of Gi exhibit wide varia- 
tions in the literature, typically Gi ~ b 2 $3 As advertised, 
Hi, ~ H s . The same situation is encountered in 2D, ex- 
cept now H s ~ H b ~ bH& L (0). The fact that the crite- 
rion for the breakdown of the Landau level-based theory 
derived from the high-field side agrees with the region 
of validity of our 3D XY-like approach derived from the 
opposite, low-field side, is another argument in favor of 
the gauge theory scenario. 

While it is the GL theory (0) that provides a realis- 
tic description of fluctuation behavior in extreme type-II 
superconductors, many numerical studies are performed 
on the 3D XY model. The ultimate low-field critical be- 
havior should be the same and computational effort is 
much reduced. It is therefore useful to discuss here the 
physical meaning of the "high" and "low" field regimes 
in the context of the frustrated 3D XY model. There is 
an immediate difference between this model and the GL 
theory ([[]) regarding the high field behavior. In the GL 
theory this regime is dominated by Landau levels and is 
characterized by strong amplitude fluctuations. In con- 
trast, in the 3D XY model, the amplitude fluctuations 
are frozen at the "microscopic" level of a single XY spin. 
As a result, there is no Landau level formation in this 
model. Instead, the high field behavior of a uniformly 
frustrated 3D XY model, as one approaches the "mean- 
field" H c %(0), is entirely determined by the pinning of 
field-induced vortices by the underlying lattice. We can 
think of this situation, to some extent, as having the LL 
structure of Eq. ^ thoroughly "mixed" by a very strong 
external periodic potential. There is, however, a relation- 
ship between the low-field critical behavior of the 3D XY 
model and the GL theory. It derives from our concept of 
"screening" of large thermally-generated vortex loops by 
the background of field-induced vorticity. In the weakly 
frustrated 3D XY model such "screening" is measured by 
a parameter k s xy, which is the XY model counterpart 
of n s in the GL theory (§): 

l f T (T,H) 1 

where measures the uniform frustration and is the 
fraction of the elementary flux quantum <j) per plaque- 
tte, while f T (T,H) is the average density per plaquette 
of vortex and antivortex segments || H piercing the xy- 
plane. Note that f T (T,H) includes all such vortex seg- 
ments, not just those connected to infinite vortex loops. 
In order for the system to be in the low-field critical 
regime of a weakly frustrated 3D XY model we need 
Ks xy»lor U<^f T (T,H). 

The actual value of H s (or Hb) in high temperature su- 
perconductors is of considerable importance. There are 
numerous estimates in the literature, based both on Eq. 
(||) and on the analysis of various experimentally mea- 
sured quantities in terms of either the GL-LLL theory 
or the so-called "3D XY scaling" (see Sec. IV). A di- 
rect estimate from Eq. (||), in a moderately anisotropic 



HTS system like optimally-doped YBCO, uses Gi « 0.01 
and Hg L (0) ps 160 Tesla, leading to H s ~ 1 - 2 Tesla. 
This estimate is subject to an irksome uncertainty, both 
intrinsic (due to our inability to theoretically determine 
H s or Hb with a precision better than within factors of 
order unity) and extrinsic (due to difficulties in extract- 
ing precise values of the GL parameters entering Eq. ([!]) , 
although the situation here is rapidly improvinga). The 
estimates of Hi, based on the fits of fluctuation thermo- 
dynamics to the GL-LLL theory are in general agreement 
with the above value of l-j2-|Teslatl and seem to give an 
upper limit Hb < 8 Tesla.e3 Similar analyses, based on 
the fits to a low-field "3D XY scaling" , generally produce 
results which seem consistent with the 3D XY-like behap. 
ior to much higher fields, 14 Tesla or even higher .EmH 
An important difference between the two approaches is 
that, within the GL-LLL theory, not only the scaling 
law but the scaling function and explicit expressions for 
thermodvMajpjc quantities are known with considerable 
accuracy.li3u'Ej In the 3D XY approach only the scaling 
law itself is known but the actual scaling function and, 
more importantly, the physics behind it are not. The 
gauge theory scenario should help remedy this situation. 

IV. CRITICAL THERMODYNAMICS AND 
^--TRANSITION 

This brings us to the second important consequence of 
description (^) , which has bearing on the nature of crit- 
ical behavior in the low- field (H <C H s ), extreme "type- 
II" limit of the gauge theory. The most significant prop- 
erty in this regime is that, for e small but finite, there is a 
true thermodynamic phase transition separating the high 
and low temperature phases of the theory, the "normal" 
and the "Meissner" state, respectively. For e = this is 
the standard H = phase transition of the Ginzburg- 
Landau theory. This phase transition is continuous and 
in the universality class of the 3D XY model. The ac- 
tual mechanism of the phase transition is directly tied 
to the expansion of thermally-generated vortex loops, as 
depicted in Fig. 2. In the ordered state below T c0 , there 
is a finite average size for such loops, A$, and configu- 
rations which contain infinite loops, "percolating" from 
one end of the system to another, do not contribute to 
the partition function in the thermodynamic limit. At 
distances much larger than A$, there is nothing to dis- 
turb the long range correlations in (<I>(0)<i>*(r)): it is not 
possible to "polarize" closed loops at such large distances 
and they behave as bound "dipoles" . This is what en- 
ables the long range phase order that characterizes the 
superconducting state. Above T c q, as more and more vor- 
tex segments are created by thermal excitation, the loops 
connect, in the sense that now there is finite contribution 
to the partition function from configurations having infi- 
nite loops, "percolating" across the system. This implies 
that A$ — > oo and it is now possible to "polarize" the sys- 
tem of loops over arbitrary large distances. Such infinite 
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vortex loops act as "free charges" and produce "metal- 
lic screening" of small external magnetic fields, resulting 
in vanishing of the helicity modulus, as discussed in the 
next section. This picture of the 3D XY phase transition 
as vortey loop "expansion" r has its origins in the works by 
Onsagepepa and FeynmarO, in the context of superfluid 
helium ,E3 but should equally well apply to high temper- 
ature superconductors (HTS) with their shortpBCS co- 
herence lengths and extremely large k (~ 100) Jl3 

As finite e (finite H in ([]])) is turned on in Eq. ([!]) we 
are facing potentially dramatic change in this picture. In 
the neutral-superfluid picture described previously, vor- 
tex loops have long range London-Biot-Savart interac- 
tions. Once e is finite, these interactions are "screened" 
by the vector potential S and, at distances much longer 
than the "penetration depth" A s oc 1/e, all the inter- 
actions are short-ranged. The simplest and best known 
example of this is just an ordinary superconductor at 
zero external field. There e is the real electrical charge e, 
while S turns into the ordinary Maxwell's vector poten- 
tial A. This charged-superfluid-|problem has been studied 
extensively, starting with Ref.EZI,_and is presently thought 
to have the following propertieseil: as already indicated 
in Sec. II, the charge e is a relevant operator in the RG 
sense, with scaling dimension equal to 1/2. This immedi- 
ately destabilizes the neutral-superfluid, zero-charge 3D 
XY critical point. There are, however, two new critical 
points, characterized by finite charge. The behavior of 
strongly type-II superconductors (k ^> 1) is determined 
by the stable critical point and describes the continu- 
ous phase transition between the normal state and the 
Meissner phase in real superconductors. El Another criti- 
cal point is tricritical and unstable in one RG direction, 
in addition to temperature. This tricritical point defines 
the transition between type-II (small charge) and type-T 
(large charge) behavior and takes place for n ~ 0.4/\/2.tll 
In a type-I superconductor, the phase transition is ex- 
pected to be discontinuous, as originally argued in Ref.LJ. 
In the type-II regime, where the transition is continuous, 
the universality class for the, pharged-superfluid appears, 
to be the "inverted 3D XY"c30, or very "close" to itO'tHl 
(see Appendix A for further details). 

What is the connection between these general proper- 
ties of charged-superfluid model and our problem? It 
stems from the gauge theory (^). This theory looks 
just like the theory for charged-superfluid, except-for the 
charge anisotropy, which should be irrelevantEjCJ. The 
underlying physics, of course, is very different. There 
is no fluctuating electrodynamic vector potential in our 
case, since we are in the k — > oo limit. Instead, our fic- 
titious vector potential S describes the long wavelength 
vorticity fluctuations in the background s vortex system 
and our "charge" e is the original magnetic field H in 
disguise (e 2 cx \/~H). Despite this difference in phys- 
ical meaning, the long distance behavior of (|J) should 
still be closely related to the electrodynamics of charged- 
superfluid. In particular, we expect two different thermo- 



dynamic phases of (|l|): the high temperature phase with 
only short range correlations in ($(0)<I>*(r)) (($) = 0) 
and the low temperature phase, in which ($(0)$*(r)) 
develops long range order (($) ^ 0). The "Meissner 
phase" (or the ^-ordered state) of fl) corresponds to the 
state of the original GL theory (P in which only 7V$ 
field-induced vortex lines cross the system from one end 
to another along H. All other vortex excitations form 
either closed thermally-generated loops of finite size or 
finite "overhang" configurations decorating field-induced 
lines as they make their way meandering from bottom to 
top of the sample. These field-induced vortex lines, or 
s vortices in reformulation (|^), have a finite line-tension 
relative to the field direction and undergo effective "diffu- 
sion" along the z-axis (this is discussed in greater detail 
in the next section). In the high temperature, "normal 
metal" phase of (f|), the $-order is destroyed by expan- 
sion of thermally-generated vortex loops and "overhangs" 
decorating s (field-induced) vortices. We now have new, 
thermally generated infinite loops "percolating" all the 
way through the system in all directions. These new infi- 
nite loops come on top of the always present backgrowpd 
of N$, s vortices. This is the nature of the $-transitionEi in 
the gauge theory ([I|) and in the reformulation (^]). The 
^-transition is the finite field version of the zero-field 
superconducting transition. a Its thermodynamics, how- 
ever, belongs to a different universality class: charged- 
superfluid ("inverted 3D XY") as opposed to neutral- 
superfluid (3D XY) at H = 0, with finite H playing the 
role of finite charge (e 2 cx VH) in the gauge theory (||). 

The ^-transition line, T$(iJ), plays a pivotal role p. 
the gauge theory scenario. Since, on general grounds,Ell 
we do not expect any true criticality associated with the 
first-order vortex lattice melting line in 3D, T$(H) is the 
only critical line in the H — T phase diagram of the origi- 
nal GL theory (Q) and controls fluctuation thermodynam- 
ics and transport at weak magnetic fields. It decides the 
issues of relevance or irrelevance of various terms that 
can be added to (Q) (point or columnar disorder, true 
electromagnetic screening with finite k, etc.) and pro- 
vides the foundation on which one can build a meaning- 
ful phenomenology of extreme type-II superconductors. 
In this respect, $(r), the new "superconducting" order 
parameter characterizing the "line liquid" state, is "more 
fundamental" than the original \l/(r). This will now be 
amply illustrated. 

To start building such phenomenology, we first need 
a reasonable estimate of T$(iJ) (Fig. 1). It starts at 
the zero-field superconducting transition T c q, where 
and \&(r) are one and the same: the H — > limit of (jl]) 
coincides with the A$ — > limit in Eq. @ and with the 
e — > limit of the gauge theory ([!]). At finite, but weak 
field, we are in the "extreme type-II" regime (n s 3> 1) of 
the gauge theory (||) and we expect that the ^-transition 
is continuous and immediately becomes "inverted" . The 
transition temperature T$(H) is gradually reduced as 
function of H (or, equivalently, e in ([|)), due primar- 
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ily to the reduction in the effective line tension of very 
large (3> £) vortex loops caused by "screening" generated 
by the "medium" of field-induced vorticity. At these low 
fields, T<s>(H) can be evaluated directly from ([|). As H 
increases, however, numerous additional terms present in 
Eq. (H) , but not included in the gauge theory ^) on the 
grounds of their RG irrelevance at long distances (3> £), 
start affecting T$ (H ) . Among such terms none are more 
important than short distance (~ £) positional correla- 
tions which eventually lead to s vortex lattice formation 
at low temperatures. 

In general, the ^-transition and the vortex lattice melt- 
ing are two completely different phase transitions, with 
two different order parameters, driven by two different 
mechanisms. One is a q — ► 0, another q ~ l/£ tran- 
sition. They are not entirely unrelated, however, since 
they arise in the same theory, ([!]) or (|^). Eor instance, 
as H -> 0, we must have T$(H) > T m (H)!B This is so 
because only in the <I>-ordered state do s votticcs in (||) 
have finite long range interactions, oc |($)| 2 .l3 Without 
such long range interactions the s. vortex system would 
remain in a liquid state as H — > O.u Similarly, in the solid 
phase, s vortices form a lattice and cannot screen large 
thermally-generated vortex loops, i.e. e becomes effec- 
tively zero even for H ^ 0. All vortex loops will then 
remain small and bound, just as they were at H = 0. 
This is discussed in more detail in Sec. VII. The prob- 
lem is that the melting transition is always first-order 
and thus, in principle, we could have T$(H) = T m (H) 
over most of the H — T phase diagram or even for all 

H. This would mean that melting is so strongly dis- 
continuous that it always "jumps" over the intermediate, 
^-ordered phase, straight into the true normal state. I 
know of no argument to rule out this possibility. 

This being said, the most likely outcome is the one de- 
picted in Fig. 1. At higher fields, as we approach the 
"type-F regime of (@) (H ~ H s (|)), the gauge theory 
suggests that ^-transition itself converts to first- order. In 
this situation, it seems justified to assume that T<j,(7J) = 
T m (H), as shown in Fig. 1. For low fields, H <C H s , 
where the melting transition becomes weakly first-order 
and (0) predicts a strong "type-II" behavior and continu- 
ous transition, it is natural to expect T^>(H) > T m (H). 
At fixed low field, as we increase the temperature in Fig. 

I, both the effective strength of the Biot-Savart inter- 
action between s vortices (Sec. VII) and their effective 
"mass" (Sec. V) decrease. As interactions and line ten- 
sion go down, a natural progression of thermodynamic 
phases follows: a solid (Abrikosov lattice), a "massive" 
liquid (<I>-ordered phase or "line liquid") and, finally, a 
"massless" fluid of unbound loops (a true normal state). 
A mean-field calculation, performed in Ref. 9, indeed 
leads to such results. I propose here a simple criterion 
which summarizes the results of such calculations and 
can be used to determine Tq,(H) and T m (H) at low fields, 
H <C H s : the s vortex lattice melts when the average size 
of thermally-generated loops, A$(T,H), reaches a frac- 
tion d m of the average distance between field-induced (s) 



vortices: A$(T,H) = d rn V2^{£\\/Z±) 1/3 £- d m ~ 0.2-0.3 
and A$(T, H) « A$(T,0) ~ (£l£||) 1/3 seems a reason- 
able estimate. Here £m| = £o±,o||KI" are tne true su- 
perconducting correlation lengths at H = 0. This results 
in the expression for the vortex lattice transition temper- 
ature in the critical region: 



t m (h) 



{d m V2^) 3 ^ 2 ^GL,± 



^ £o_L ) 3 / 2 fe 3/4 



(10) 



where the temperature is measured relative to the true 
zero-field superconducting transition, t = (T/T o) — 1, 



h is defined below Eq. 



and 



was set to 2/3. 



The ratio £o±/£gl,± should be ~ 1. As argued above, 
we expect d m ~ 0.2 — 0.3. The ^-transition, on the 
other hand, takes place when the size of thermally- 
generated loops, at finite H, reaches the sample dimen- 
sions, A<j,(T, H) — > oo. This should take place along the 
line where the average loop size for H = 0, A$(T, 0), be- 
comes of the order of average distance between s vortices, 
i.e. A$(T,0) = (Z$\/27r(£||/tL) 1/3 ^ with d$ - 1. This 
determines the vortex loop "expansion" line, or Tq>{H): 
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0_L \3/2 



7= ( J ^ ± ~) 



(11) 



with (2$ ~ lB Obviously, T^{H) > T m (H) since d m < 
d$. Eqs. ( p~o| , p~T|) are valid only in the limit of low 
fields, H < H s (g). At higher fields, H ~ H s , T$(H) 
and T m (H) merge together and both vortex loop "ex- 
pansion" and vortex lattice melting occur simultaneously 
when A$ (T, H) reaches ~ ^/2tt£ from within the solid 
phase (Fig. 1). The above expressions ( p^|Jll| ), with d m 
and d$ serving as numerical parameters, can be viewed 
as a "Lindemann criterion" for vortex loops and should 
provide reasonable estimates of X$ (H) and T m (H) . The 
reader should bear in mind that the above expressions for 
T m (H) and T$(H) are only estimates, however reason- 
able. Only the first-principle calculation, starting from 
(|l|) or (|J), and computing accurately free energies of dif- 
ferent phases (solid, ^-ordered liquid, and the true nor- 
mal state) can settle this question within the theoretical 
framework proposed in this paper. Unfortunately, such 
calculation is not feasible at the present level of analyt- 
ical sophistication. Additional careful experiments on 
HTS and numerical simulations of 3D XY and related 
models are far more promising in this regard. 

I now proceed to further investigate the phase diagram 
represented by Fig. 1 and Eqs. ( |To| , pT| ). As the field is 
turned on in Eq. (Q) we can immediately write down the 
scaling expression for the dimcnsionlcss singular part of 
the free energy, /, associated with critical fluctuations: 



f=\A 



H 



H k t 



-) 



(12) 



where t = (T/T c q) — 1 and Hk depends on material pa- 
rameters. This expression is completely general and as 
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such conveys little information. It is based only on the 
existence of the zero- field critical point. The same expres- 
sion can be written for spin systems or any other system 
exhibiting a critical point which is then perturbed by 
a generalized "field" .E3 In our case, the H = critical 
point is in the 3D XY universality class and we should 
have a — a xy . Furthermore, based on dimensional anal- 
ysis r^pd general physical arguments, it was proposed in 
Ref.O thai—A = 2v xy . This result holds to two- loop or- 
der in RGc3 and is likely exact, as emphasized in Sec. II. 
4>±{ h ^a ) is a universal function of its argument inside 
the 3D XY critical region of the GL theory. At present, 
its form is not known. Note that Eq. (|12|), while com- 
pletely general, is written in the form which implicitly 
suggests that the finite-field critical behavior is governed 
by the zero-field critical point, as is frequently the case 
in spin systems. 

In the gauge theory scenario, the situation is different 
and we can be more specific. First, as already empha- 
sized, within this scenario the critical fluctuations are 
governed by a critical line and not a critical pointn This 
means that we immediately learn something about the 
function <f>± defined in Eq. ( |12|): 6 + is non-analytic along 
the ^-transition line, T^{H) (|ll|). This line singularity 
should be explicitly incorporated into the expression for 
the free energy. To devise such a new scaling function, 
based on the gauge theory scenario, we start by observ- 
ing that we can eliminate the "trivial" part of the charge 
anisotropy from (^) by rescaling all lengths and fictitious 
vector potential S with an appropriate superconducting 
correlation length, £x,||> in the way that makes the 
dependent part of (|4|) isotropic. This rescaling procedure 
is a variation on the familiar rescaling of anisotropy at the 
H = transition. After the rescaling, the ^-dependent 
part of (Q) describes an isotropic superconductor with a 
correlation length £ = (£x£||) > while the coupling con- 
stants in the last two terms become: 



K 



-Ml 



^i,ii = c±,n(Cii/a) 1/3 ^ 



(13) 



The following quantities appear in Eq. (^3|): = 

£ 0U \\t\- v and £ = (£l£||) 1/3 = &\t\- v are the true di- 
verging superconducting correlation lengths at T c(h de- 
fined by the eigenvalues of the helicity modulus tensor 
(see the Appendix B). Accordingly, T = is tne <rMe 

anisotropy ratio at the H = critical point (T c q) and not 
the GL anisotropy, Tgl — £gl,x/£gl,\\- It now becomes 
clear why K± n have been defined in Eq. (|J) with the 
anisotropy T explicitly factored out: new rescaled cou- 
pling constants can simply be written as: K' ± = c± t \\T£, 
where 



pl/3 



corresponds to H — T 2 ^ 3 H 



(14) 



H is just the rescaled magnetic field appearing in the 
original GL theory (|l|) after the anisotropy at the H = 
critical point has been rescaled out. Consequently, cx,|| 



describe the fundamental anisotropy of the gauge theory 
([|), which is inherent to the H 7^ problem a»d-4s not 
associated with a "trivial" anisotropy at T c o£3'E3 The 
corresponding fictitious "charges" associated with K' ± 
are: 



e 2 



(T,H) = 



1 



c±,\\(T,H)e 



(15) 



The product of the above rescaling procedure is a ficti- 
tious anisotropic electrodynamics with two "charges" , ej_ 
and en. As discussed in Sec. II, the charge is a relevant 
operator at the H = (e = 0) critical point, with scaling 
dimension equal to 1/2. We thus define two dimension- 
less scaling variables: 



e 2 



■JL.II' 



t>o cx.ii \t\ v 



(16) 



Note that dimensionless ratio £/£ is common to both 
charges and (non-trivial) anisotropy is stored in c± and 
cm. c± and cm, however, are also functions of £/£ only. 
Therefore, there is only a single relevant scaling variable, 
the dimensionless charge: 



% 2 



/2^ r i/ 3 tf 1/2 



(17) 



which is precisely the original scaling variableEZI of the GL 
theory (^, since £/£ = £x/l- The functions cx,||(£/^) are 
discussed further in the next section and Appendix B. 

We are now in position to write down the scaling func- 
tion for the free energy within the gauge theory sce- 
nario, with the non-analytic part associated with the $- 
transition explicitly factored out: 



fs = \t-U(h)\ 



f t - t^(h) \ 
\ Mh)\ J 



(18) 



defined in Eq. (|TT[) , also follows from 
is therefore a universal number of the 



Note that t$(h), 
q 2 = \/27rd$. c?$ 
GL theory (|j), as is d m £2 fi^' (x) is a universal and reg- 
ular function of its argument. The subscripts ± refer to 
x > (x < 0), while the superscripts L and S indicate the 
"vortex liquid" and "vortex solid" branches of f2, respec- 
tively. For example, below T m (H), we should use il^_(x). 
In writing down ( pj| ) I have assumed that the correlation 
length exponent of the gauge theory vqt ~ v xy ~ 2/3 



and that the hyperscaling relation holds, resulting in 

How do we evaluate the crossover function, fl(x)7 
I alert the reader to the following important point: 
the gauge theory scenario explored in this paper allows 
one to, in principle, determine all the branches of the 
crossover function f2(x) by using a combination of pertur- 
bation theory and RG techniques. Such analytic calcula- 
tion is extremely laborious and far beyond the scope of 
this paper. Well informed reader will immediately realize 
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that many aspects of this calculation are computation- 
ally extremely demanding, and actually have not been 
accomplished in the published literature even for the or- 
dinary H = situation. Indeed, the technical difficulties 
involved are of the same general nature. This, however, 
does not detract from the main message of this section: 
the underlying physical picture of the gauge theory sce- 
nario provides a systematic, conceptually straightforward 
way to compute the H ^ 3D XY critical thermodynam- 
ics at the same level of analytical accuracy as is presently 
feasible for the H = case. 

Faced with such odds, I assume, for the purposes of this 
paper, that Q(x) in Eq. |[|) is some unknown universal 
crossover function, to be determined either from numeri- 
cal simulationsEa or directly from experiments. With the 
free energy thus specified, we can proceed to evaluate the 
singular part of all thermodynamic functions, simply by 
taking requisite derivatives .c3 



V. HELICITY MODULUS, LINE "DIFFUSION", 
TOPOLOGICAL WINDINGS, AND PHYSICAL 
NATURE OF $-ORDER 

I now turn to physical properties which allow more di- 
rect look at actual configurations of loops and lines that 
characterize the state of a superconductor above and be- 
low T$(iJ). A useful measure of a degree of supercon- 
ducting order is a q-dependent helicity-modulus tensor, 
T(q), whose components are defined asL3 



V(q) = V 



S 2 F 



(5a 1/ (q)(5a M (-q) 



(19) 



where fj,, v = x, y, z, V is the total volume, F is the 
free energy of the GL theory ([!]) and a(r) is a small 
(infinitesimal) vector potential added to the external A. 
The uniform component of the associated magnetic field, 
h(r) = V x a, vanishes. The above second derivative is 
evaluated in the a — * limit. 

T(q) measures the ability of a system to "screen" 
out tiny external fields. In the superconducting phase 
lim q _»o Y(q) is finite and the system is said to exhibit a 
differential Meissner effect. In a normal phase, T(q) ~ q 2 
and vanishes in the q — > limit. Within our "helium 
model" , the way T is reduced to zero in the long wave- 
length limit is through proliferation of infinite vortex 
loops/lines which go all the way across the system and 
can act as "free charges" , screening a weak external per- 
turbation. In this intuitive sense, we can think of a nor- 
mal state as a vortex "metal" , while the superconducting 
state is a vortex "dielectric" , with only vortex loops of 
finite size present as thermal excitations. 

To compute the helicity modulus of our original GL 
theory one adds an infinitesimal a„ to A M in Eq. (|l|). If 
we now go to the reformulation (0) and finally, through 
the coarse-graining procedure of Appendix A, end up 
with our fictitious gauge theory, a p appears as a small 



addition to the "vector potential" S M in the second (gra- 
dient) term of Eq. (0). This implies that the long wave- 
length (q <C lit) form of the helicity modulus of the 
gauge theory (Q) coincides with that of the original GL 
theory ([l9|). Using the gauge theory ffl) and ignoring the 
anisotropy, we obtain that below T$ (ff) (see Appendix B 
for details): 



T(q) - Kq 2 



K 2 q 4 
T 



(S(q) ■ S(-q)) = Kq 2 + 0{q 



(20) 



In the ^-ordered state our fictitious gauge field is "mas- 
sive", i.e. exhibits a Meissner effect, and lim g _,.o(S(q) • 
S(— q)) ex K^)! -2 goes to a finite value. The simple 
physics behind this is that thermally-generated vortex 
loops in $(r) have average size that is finite and do not 
contribute at all to T(q) in the q — > limit. Further- 
more, Eq. (^0|) tells us that, if -Kj. n are finite, the <f>- 
ordered phase is not a superconductor and has a finite 
superconducting correlation length, both perpendicular 
and parallel to the external field (see Appendix B for 
details) , 



£/£||~Aji/T. 



(21) 



This result is easily understood: with H finite, there now 
must be iV$ field-induced vortex lines moving about in 
the sample. Unless these (s) vortices are pinned down, as 
it happens in the vortex-solid phase, they will be avail- 
able to "screen" weak (infinitesimal) external fields and 
the system is always a vortex "metal" with the "screen- 
ing length" K/T. In particular, according to our as- 
sumption ii), the system of field- induced vortex lines also 
contains windings in the xy plane and such "screening" , 
although anisotropic, is finite in all directions. Below 
T$(i7), where s vortices are exclusively responsible for 
the vanishing lim g ^o T Al (q), cjii(<?q) determine the _L 
and || "screening lengths" in the s vortex system, in units 
of magnetic length t (Eq. ( pi] ) in Appendix B). 

Above T$ (H ) the situation changes and thermally gen- 
erated vortex loops "expand" across the system. Obvi- 
ously, the helicity modulus still vanishes, but now there 
is an abrupt drop in the coefficient of the q 2 term: 



T(q) = (K — G 



K 2 , 



K(l - C 



H 



(22) 



where t(T,H) = (T - T$(H))/T cQ and C and G are 
constants, vqt is the thermodynamic exponent of the 
Meissner transition in our fictitious electrodynamics (Q), 
and vgt ~ v-xy ~ 2/3, as argued in Sec. IV. The sec- 
ond term in the above equation arises from lim 9 _,o(S(q) ■ 
S(-q)) = G/^q 2 in Eq. ©, right above 
This implies that the original superconducting correla- 
tions, measured by (\E f (0)\E f *(r)}, remain finite in all di- 
rections on both sides of T$(H), but there is a non- 
analytic drop in the superconducting correlation length 
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at the ^-transition, as thermally generated loops prolif- 
erate through the system and additional infinite vortices 
become "free charges" and available to screen. The new 
order parameter, ^(r), however, does attain a true long 
range order below Tq>(H), i.e. — > oo as T — ► T<j,(H) 
from above. It is unfortunately rather difficult to mea- 
sure the ^-correlations directly, by probing some suit- 
ably defined "helicity modulus" associated with the <E>- 
order. This would require defining quantities which are 
configuration-dependent and highly non-local, a rather 
time-consuming proposition in a typical numerical simu- 
lation of a 3D XY or related model. 

Still, the situation is far from hopeless. We can devise 
another set of criteria that are relatively easy to imple- 
ment in numerical simulations and yet allow for a rather 
intimate look at the <J>-order and what precisely takes 
place as we cross the ^-transition line. Below T$(iJ), 
thermally generated vortex loops are bound and field- 
induced vortices execute an effective "diffusive" motion 
along the field direction. An average transverse displace- 
ment of a single field-induced vortex line from the point 
where it starts at z = to its ending point at z = L z 
goes as 



thermally-generated vortices proliferate across the sam- 
ple in all directions, D s — > oo (m s — > 0) and we get 
hyperdiffusion: 



P = 



(23) 



where D s is the effective "diffusion" constant. This is 
shown in Fig. 3. The cutting/reconnecting of vortex lines 
does not affect this diffusion process except by renormaL 
izing D s , as long as we are in the <&-ordered phase.Ej 
For example, in the 3D XY model, where the identifi- 
cation of an individual field-induced line is not unique, 
we should simply look at the distribution of distinguish- 
able vortex paths, obtained by randomly resolving all the 
crossings, and average over all distinct configurations. 
Such distribution will be "diffusive" , with the average 
rms displacement given by Eq. (23). We can use this 
effective "diffusion" constant D s (23) to define an effec- 
tive "mass." in the elegant non-relativistic boson analogy 
of NelsonB. Note that the worldlines of such flux-bosons 
do not correspond to (s) vortex lines in individual con- 
figurations of an extreme type-II superconductor. This 
is clear since non-relativistic bosons describe strictly di- 
rected lines, i.e. contain no "overhang" configurations as 
they advance from bottom to top of the system along z- 
axis (the "time" axis in the boson analogy) . Such "over- 
hang" configurations, plus numerous vortex loop exci- 
tations floating around, describe worldlines of "particle- 
antiparticle" creation processes and cannot be accommo- 
dated within the non-relativistic quantum boson anal- 
ogy. Still, as long as we are in the <I>-ordered phase, it is 
only the iV$ field-induced vortex lines that go all the way 
across the system. We can then define an effective system 
of N$, flux-bosons in the boson analogy, with suitably ad- 
justed bare mass m s and effective interactions, so that its 
long distance (3> £ and 3> A$(T, H)) behavior faithfully 
emulates an extreme type-II superconductor (Appendix 
A). Above T<s>(H), as infinite tangles of field-induced and 



LZ 



P 



1 



(24) 



This hyperdiffusion arises through processes depicted in 
Fig. 3, where a vortex line winding along the field direc- 
tion simultaneously winds all the way in the xy-plane by 
connecting itself to thermally-generated tangles, which 
are naturally present in the ^-disordered phase. In the 
3D XY model, this implies that the distribution of trans- 
verse displacements of individual field-induced vortex 
paths is no longer "diffusive" and has rms displacement 
limited only by the system size ([24|); more precisely, the 
distribution of acquires a power- law tail above T$ (H) . 
Such windings in the xy-plane are plainly in evidencejn 
the recent numerical simulations of Nguyen and SudboO, 
right above their melting line. With such additional xy 
windings present with a finite weight in the partition 
function, the "effective mass" m s of non-relativistic flux- 
bosons vanishes since the vortex line tension relative to 
the field direction has gone to zero. Above T$(H), infi- 
nite vortex paths of length ~ L 2 (assuming L± = Lu = L 
and a simple random walk) crossing the system in all di- 
rections contain a finite fraction of all vortex segments: 
these are the "massless excitations". In this respect, the 
^-transition corresponds to the restoration of "relativis- 
tic invariance" in a dual system of quantum particles 
whose worldlines are our original vortex loops and lines. 
To wit, the ground state of spjh quantum system, con- 
taining only "vortex matter"l£j below T$(ff), explodes 
with "vortex matter" , "vortex antimatter" and "vortex 
tachyons"Eil (Fig. 3), as the "vacuum" becomes unsta- 
ble at 7$ (H) to spontaneous creation of "particles" (Ap- 
pendix A). 

The above connection between the "$-order" and the 
effective line tension of field-induced (s) vortex lines re- 
veals directly the physical content of the gauge theory (^) 
and permits us to construct a purely geometrical picture 
of the ^-transition. To do so, consider once again the nor- 
mal state of an extreme type-II superconductor or a 3D 
XY model at H = 0. Above T c q, we can find vortex paths 
that go all the way from one end of a sample to another, 
in any direction (Fig. 2). If we work with periodic bound- 
ary conditions in all directions, this statement means that 
we have some windings along x-, y-, and z-axes. To 
understand what is precisely meant by such windings, 
we wrap our system on a three-dimensional generalized 
torus, embedded in a four-dimensional space This is just 
a geometrical way of representing periodic boundary con- 
ditions. Now, the number of windings along, say, x-axis, 
Af x , is the total number of continuous vortex paths in 
the whole system that wind all around the torus in the 
x-direction, irrespective of their orientation, i.e., whether 
they are "vortex" or "antivortex" paths relative to the 
x-axis. Such paths are topologically distinct from fi- 
nite closed vortex loops: the latter can be continuously 
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shrunk to a point while the former can not. Af x is dif- 
ferent from the winding number, W x : W x also counts all 
the windings along x-axis but with a single "vortex" path 
contributing +1 while an "antivortex" path counts as —1. 
In the widely recognized vocabulary of the 2D XY model, 
Af x would correspond to the total number of vortices plus 
antivortices in the yz-plane, while W x would be the total 
number of vortices minus the total number of antivor- 
tices. Back in 3D, in the superconducting state below 
T c o, all thermally-generated vortices come in the form of 
finite closed loops and both Af x = and W x = 0. Above 
T c o, W x must remain equal to zero due to the "vortex 
neutrality" of the GL theory (0) or the 3D XY model, 
but Af x is now finite and Af x oc where L is the linear 

size of the system (we are assuming L±_ = Lu = L) and 
u is the "anomalous dimension" of such infinite paths. 
This implies that there is a finite fraction of all vortex 
segments contributing to such windings along x in the yz- 
plane. The same holds for the winding number and the 
total number of windings along y- and z-axes, W y ( z ) — 
and Afy( z } cx L 1+u . This can be summarized as: 

Af x ,y,z = for T < T cQ , 

Af x , y>z = 2nl y , z L 1+u for T > T c0 , (25) 

where nf\ is the "density" in the yz-plane of thermally- 
generated infinite vortex-antivortex winding paths 
traversing the system along the x-axis, and so on. In the 
isotropic case = Uy = nj. Of course, the presence of 
such windings in all directions is the reason why the ma- 
terial is not in the superconducting state above T c q: these 
infinite vortex paths can now move to "screen" weak ex- 
ternal fields, driving the helicity modulus to zero in the 
long wavelength limit and producing finite dissipation. 

As we turn on a finite field in ([!]), we still have 
W x ( y ) = but W z = 7V$ and consequently Af z must be at 
least N$ in every configuration of the systemjj Imagine 
now how the state of the system evolves along a small 
circle in the H — T phase diagram (Fig. 1), surrounding 
T c q. Our circular path starts at H = and at some tem- 
perature T slightly above T c q, and evolves in the counter 
clockwise direction toward its end point at H = and 
some temperature slightly below T c0 . Initially, we are 
very close to the H — normal state and it is safe to 
assume that N^L < Af z (T,H = 0)L 2 ~ U at some tem- 
perature T slightly above T c q. In this case it is natural 
to expect that, in a finite field: 

Af x . v = 2nl y L 1+Ul - , Af z = n^L 2 + 2n T z L 1+u w , 

(26) 

where n$ = N$/L 2 = l/2ir£ 2 is the density of field- 
induced (s) vortex lines and n x Z (T,H) are "densities" 
of thermally-generated vortex-antivortex windings. Note 
that now 2nf\ ^ n$ + 2n^ and nf\ ^ rff, even in 
the isotropic case (although the gauge theory scenario 
strongly suggests u± = ity). There is a finite "density" of 



thermally generated infinite vortex (or antivortex) wind- 
ing paths in any direction. Eq. ( |26| ) describes how the 
normal state of the system (^5|) has changed after the 
application of a finite, but weak external field. 

On the opposite side of our imaginary cirle, near its end 
point at H = and T < T c q, the situation is completely 
different. Now, the zero-field state is a superconductor 
and Af x .y tZ — 0. Any finite field, no matter how small, 
has a drastic effect. For very low fields, it is natural 
to expect that there are no thermally-generated infinite 
vortex loops and only those windings associated with the 
field-induced vortex lines are present in the system: 

Af x . v =0 , Af z = A$ = n$L 2 . (27) 

This is just the (s) vortex lattice state in Fig. 1. 
Note that in these general geometrical terms there is 
no difference between the (s) vortex lattice state and 
the 'ianisotropic superconducting liquid" of Feigel'man 
et alE3. Due to the absence of windings in the xy-plane 
both have superconducting response along the field di- 
rection, i.e. the T 22 (q) component is finite in the q — > 
limit. Also, in both cases, T^x^) (q) vanish as q — > 0. 
The only difference is that it takes arbitrary weak pin- 
ning to restore superconductivity in all directions in the 
vortex lattice state. On this basis, I have assumed in the 
phase diagram of Fig. 1 that such an "anisotropic super- 
conducting liquid" phase is preempted by the first-order 
transition at T m (H)£3 

In our proposed phase diagram depicted in Fig. 1, 
the intermediate ^-ordered phase is inserted between the 
true normal state (^6|) and the vortex lattice phase (|27|). 
What is its nature in simple geometrical terms used to 
describe the other two phases in Eqs. ( p6| ) and (f27|)? In 
the <&-ordered state all thermally-generated vortex loops 
are bound and only iV$ field-induced (s) vortex lines go 
from one end of the sample to another, along H, result- 
ing in precisely N§, windings along the z-axis. However, 
these (s) vortex lines are in a liquid state, characterized 
by some finite effective line tension, T, and they "diffuse" 
in the xy-plane while winding along H. This translates 
into a non-vanishing total number of windings in the xy- 
plane: 

Af x . v = I x , y L 2p , Af z =N< s> =n$L 2 . (28) 

We expect 1 > p = 1/2.S I x<y are some finite quantities 
having dimension of (length) ~ 2p and we use Eq. (E8|) as 
their definition. Eq. ( p8| ) describes the <I>-ordered state 
using a simple geometrical language of this section. The 
fact that Af x ,y 7^ has nothing to do with the thermal 
"expansion" of vortex loops; these are still all of finite 
size. Rather, it is due to the lateral "diffusion" of field- 
induced lines, as they wind along the z-axis. The field- 
induced vortices tend to form infinite "clusters" , which 
manage to wind a finite number of times, Af xy , along 
x(y) by winding infinitely many times, Af z ~ L, along 
the field direction (I now set p = 1/2). For such infinite 
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clusters, we can define a long distance, effective "diffu- 
sion" constant, 2?, and the associated line tension, T, 
both relative to the field direction: 

J\fc £ 

Defined in this fashion, 2? and T are truly global quanti- 
ties, detached from complicated individual configurations 
of interacting vortex loops and lines and dependent only 
on thermodynamic state of the system as a whole. In the 
^-ordered state both are finite and have a well defined 
average. In the true normal state, however, the infinite 
clusters appear for which 2? — ► oo and T — ► 0. Infi- 
nite vortex paths inside such clusters manage to wind in 
the x(y) direction after only a finite number of windings 
along the field. In this sense, T can serve as a probe 
for the presence/absence of the $-order. Above T$(H), 
such clusters with T = contain a finite fraction of all 
vortex segments. Again, in the <&-ordered state of the 
3D XY model, where there is no unique identification of 
individual vortex paths due to their crossings, every con- 
figuration that contributes to the thermodynamic limit 
has precisely 7V$ vortex lines going from bottom to top in 
every distiguishable assignment of such paths. All other 
vortex paths form finite closed loops. The distribution 
of the total number of windings, obtained by randomly 
resolving all the crossings, still satisfies (pq). 

Once the system makes this phase transition to the 
normal state ( p6[ ) and the $ order is lost, there is 
only a smooth evolution. The "densities" of thermally- 
generated windings n^ yz , as well as the way these wind- 
ings are realized in individual configurations of vortex 
loops and lines, can change considerably, depending on 
where we are in the H — T phase diagram relative to 
T$(i2), but we do not expect to cross any additional 
phase boundaries. 

The geometrical picture presented here, based only on 
global topological properties of loops and lines, gives a 
clear insight into two forms of "superconducting" order, 
described by two different order parameters, W(r) and 
$(r). The familiar superconducting order, measured by 
^(r), reflects the system's ability to expel tiny exter- 
nal fields and is manifested by a finite helicity modu- 
lus and absence of linear dissipation. This leads to the 
well-known spectacular experimental consequences and 
is naturally of great practical importance. A more sub- 
tle form of "superconducting" order, associated with the 
new order parameter ^(r) introduced in RefB, describes 
the presence of finite line tension at all lcngthscales and 
is manifested by suppression of large thermally-generated 
vortex loops in the partition function. At H — 0, or in 
an infinitesimal field with only a single field-induced line, 
VP and $ are equivalent. The helicity modulus and the 
effective line tension vanish simultaneously at a single 
superconducting transition, T = T c q. In a finite field, 
with a finite density of field-induced lines, the situation 
is different: the "superconducting" transition can now 



be viewed as split in two branches, T m (H) and T$(H). 
At T m (H) the standard superconducting order is lost as 
the vortex lattice melts into a liquid of field-induced (s) 
vortex lines. Even though the total number of windings 
along H is still locked at iV$, just like in the vortex lattice 
state, the effective "diffusion" of (s) vortex lines leads to 
windings in the xy-plane, N x . y oc L 2p . This amount of 
winding suffices to cause vanishing of the helicity mod- 
ulus (|2C]) and drives *(r) to zero (i.e., (*(0)^*(r)) is 
short-ranged). In contrast, the long distance line tension 
T oc 1/2? ( |29[ ) is still finite below T$(i2), same as in 
the Meissner state of the H = superconductor. This 
state differs from the normal state by the presence of 
long range $-order, measured by ($(0)$*(r)). It can 
be considered isomorphic, in the long wavelength limit, 
to a liquid (or solid, as shown in Fig. 1) of iV$ field- 
induced vortex lines with some finite effective line ten- 
sion relative to the field direction, and interacting via a 
long range, London-Biot-Savart-type interactions, whose 
overall strength is set by ~ |(<I>(r))| 2 .EJ Only at some 
higher temperature T$(iJ), does T vanish and infinite 
vortex loops proliferate across the system. This signals 
the destruction of the <3>-order as the system finally makes 
transition to the true normal state. Above T$(i2), both 
^ and $-order are absent. Consequently, the true nor- 
mal state of the GL theory in a finite magnetic field (|l|) 
should not be jdewed as a "line-liquid" in the commonly 
accepted senseEJ. The general geometrical arguments of 
this section preclude such identification and point to a 
direct connection between the presence (absence) of the 
<I>-order and our ability (inability) to emulate the long 
wavelength behavior of the system ([!]) in terms of a con- 
ventional "line-liquid" (or a "line-solid" ) of field-induced 
vortices. In the language of the non-relativistic "boson 
analogy", the mass of the bosons vanishes at T^(H) and 
such analogy breaks down in the true normal state. 

The above criteria should allow one to clearly distin- 
guish the high and the low temperature states of the sys- 
tem in numerical simulations and to establish whether 
they are separated by a true thermodynamic phase tran- 
sition or a sharp crossover. Such procedure should be su- 
perior to measurements of specific heat which possesses 
at most nearly-logarithmic singularity and is severely lim- 
ited by the finite size effects. 



VI. FLUCTUATION CONDUCTIVITY 

Electrical conductivity (or resistivity), is a quantity 
easily measured experimentally. Unfortunately, fluctua- 
tion conductivity, while in principle a very useful probe 
of a degree of superconducting order, is not purely ther- 
modynamic quantity and cannot be evaluated from the 
GL theory (|lj) unless additional assumptions are made 
concerning the time-dependence of the fluctuating super- 
conducting order parameter. We adopt here a frequently 
used and empirically successful assumption of dynamical 
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scalingj23 which connects the decay of spatial correlations 
with that of time correlations: r sc ~ £* c , where r sc is the 
relaxation time associated with the superconducting or- 
der parameter and z is the dynamical critical exponent. 
At zero field, this assumption leads to the expression for 
the fluctuation conductivity: 

o- « €lc 2 ~ D ~ t u *y( D - 2 - zS > , (30) 

with z ~ 1.5 from numerical simulations.il For simplic- 
ity, we have suppressed anisotropy in the above expres- 
sion. When a finite external field is turned on, and we 
are inside the critical region near T c q, where £ sc is very 
long, we can still write: 

aoc % C (B ^0)] 2+2 ~ D . (31) 

Strictly speaking, the dynamical exponent could be dif- 
ferent in the finite field case but I will ignore such possi- 
bility. I also concentrate on dissipative transport and do 
not consider Hall conductivity. Right below T$(i7), the 
gauge theory (Q) suggests that £ sc is finite in all direc- 
tions and £ sc oc K/T ~ c£, as discussed in the Appendix 
B. Above T<j,(H), the screening length A, defined in the 
Appendix B by Eq. (|||), drops abruptly ( f2l| ) and we 
assume that the superconducting correlation length £ sc 
does the same, 

c£ t vgt 
Wxc£(l-G—) = c£(l-C^=) . (32) 

This assumption seems perfectly justified on physical 
grounds, since the drop in the screening length arises 
through appearance above T$ (H) of thermally-generated 
infinite vortex loops which lead to additional screening. 
It is natural to expect that these same loops suddenly 
increase dissipation and produce a non-analytic drop in 
the conductivity. 

After restoring the anisotropy, the relative change in 
the conductivity, Scr^, as one crosses over from the i>- 
ordered to the true normal state is: 

6a, = ^< = = (z + 2 - D)CJ^L , (33) 

where p = (_L, ||), cr^.<.> is the fluctuation conductiv- 
ity below (above) T$ (H) , are constants depending on 
material parameters, and vqt ~ v xy ~ 2/3. The vortex 
loop "expansion" that takes place at T$ leads to a non- 
analytic increase in dissipation and a corresponding drop 
in conductivity. 



VII. VORTEX LATTICE MELTING IN THE 
CRITICAL REGION 

A theory of the vortex lattice melting in the critical 
region is an elaborate subject and its detailed discus- 
sion will be presented elsewhere. There are, however, 



several important consequences of the present gauge the- 
ory scenario that concern the very nature of the melt- 
ing transition. We therefore discuss here the "minimal" 
set of requirements that should be satisfied by any the- 
ory of melting consistent with the gauge theory scenario. 
We should first observe that the reformulation obvi- 
ously has a (s) vortex lattice as its ground state at low 
temperatures, just as the original GL theory (Q). The 
gauge theory (J|), however, does not, for the simple rea- 
son: up to this point we were interested in the long wave- 
length, q <C l/£, behavior. The ^-transition, for exam- 
ple, is clearly the q —> transition. On this basis, we 
have dropped a large number of terms from Eq. (|4|), by 
arguing that they are irrelevant at very long distances. 
The melting transition, in contrast, is a, finite g-transition 
(q r-j l/£), and it requires additional terms and modifica- 
tions to the coarse-graining procedure applied on the way 
from (||) to (|]). For instance, higher powers of S, partic- 
ularly the odd ones, reflecting the up-down asymmetry 
along H manifest in Eqs. (|l]) and (0), where unimpor- 
tant in the long wavelength limit but are essential for the 
transition to a non-uniform (s) vortex lattice state. 

This being so, the ^-transition casts a long shadow 
on the melting transition. This is clear from (Q) and, 
by inference, from (^). In the <I>-ordered phase, there 
is an effective long range interaction between (s) vortex 
lines which, in. the q -C 1/1 limit, takes the London-Biot- 
Savart form £3 To see this, note that, in the <I>-ordered 
"Meissner" phase, our fictitious "photon" is "massive", 
i.e., the second (gradient) term in Eq. (Q) becomes: 

Tx |($)| 2 Si+7||l($)| 2 Sl , (34) 

where ($) is the order parameter associated with the <&- 
order. As seen in Sec. V, we expect that, as long as 
($) is finite, we can write some effective description of 
such a state in terms of a "lineJiquid" of field-induced 
vortices, in the sense of NelsonE Long distance physics 
can be described through fluctuations in the density and 
"currents" of vortex lines, Sp(r) and j(r) — {j x ,j y ), 
respectively.H The connection with our fictitious "gauge" 
potential S(r) (||) is (Appendix A): 

(V X S) || ->■ 2n5p , (VxS) ± ^2ttJ . (35) 

If we now reexpress (|34j) in terms of Sp and j we get pre- 
cisely the long distance part of the London-Biot-Savart 
interaction between field-induced vortex lines: 

4AlW| 2 V ' P(q)M ' q 2 )+ ^ q) - jVq) , (36) 

q ?± + q \\ 

where the continuity condition, q- j(q) = —q\\Sp(q), is as- 
sumed. The anisotropy, suppressed in the above expres-. 
sion for simplicity, can be straightforwardly restored.o 
This expression ( pq ) for the effective interaction at long 
distances is just what is obtained in the mean-field based 
approach,Q but with one crucial difference. The overall 
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strength of the interaction is not given by the mean-field 
amplitude squared of the superconducting order param- 
eter, but by |($)| 2 cx n$, where n$ is the $-superfluid 
density, whose physical meaning is apparent from refor- 
mulation (||) and gauge theory (Eh. 

The immediate consequence of (^6|) is that the melting 
line, T m (H), goes into T c o, the true zero-field supercon- 
ducting transition, as H — > 0. This result is strongly 
suggested by all available numerical simulations on the 
3D XY model and arises naturally in the gauge theory 
scenario. For all its apparent simplicity, this result is not 
trivial: the mean-field based theories .of melting includ- 
ing only field-induced London vorticesQ naturally lead to 
T m {H) j=^_T c , the mean-field transition temperature, as 
H — > 0E3E3 Furthermore, the exponent v has its mean- 
field value t'mf = 1/2 and is not equal to v xy ~ 2/3. 
Therefore, the thermodynamics of the melting transition 
resulting from such theories cannot satisfy the 3D XY 
scaling properties of Sec. IV. This is a direct consequence 
of ignoring those very degrees of freedom (vortex loops) 
which are primarily responsible for moving the true zero 
field superconducting transition temperature from T c to 
T c o and changing v m { to v xy in the first place. This is a 
serious flaw and must be rectified in a proper theory of 
vortex lattice melting in the critical region. 

As a first step, we attempt to remedy the situation 
by simply replacing, by hand, the mean-field amplitude 
squared with the true superfluid density at zero field. 
This amounts to installing \(&)h=o\ 2 instead of K^)) 2 in 
Eq. (|36|). With the interaction fixed in this fashion, we 
can then proceed to analyze the same, finite field model 
of Refill, still including only the field-induced vortices. 
This is equivalent to having T 9 (H) = T^(H = 0) = T c0 , 
as represented by the dashed vertical line in Fig. 1. 
This procedure, arbitrary as it is within the framework 
of mean-field based theories, seems to get us out of the 
above difficulties, since now evidently T m (H) — > T c q for 
H — > and the exponent v takes its true H — value, 
v xy ~ 2/3. 

However, things are not that simple: to understand 
why note that the above remedial procedure is in fact 
exact, but only for a single field-induced line. For a 
finite density of lines, the physical state of thermally- 
generated vortex loops and other fluctuations described 
by <&(!■), which controls the effective interaction between 
field-induced vortex lines through \{<&)\ 2 in Eq. (|36"|), is 
itself strongly affected by interactions with those same 
field-induced lines. The effective coupling of these two 
interpenetrating systems, vortex loops and lines in refor- 
mulation (Q), must be solved self- consistently at finite H: 
this is precisely what is accomplished in the gauge theory 
scenario (Q) for the long wavelength (q <C l/l) behavior. 
Clearly, a "minimal" theory of vortex lattice melting in 
the critical region must involve both the positional order 
parameter of the vortex lattice pa (or the original "J) and 
the new "superconducting" order parameter <!>. The cou- 
pled equations governing the (T, H) dependence of these 
two order parameters must be solved simultaneously and 



self-consistently near T m (H). 

An important physical feature is expected to emerge 
from such a solution: the formation of vortex lattice is a 
phase transition involving simultaneous ordering of both 
field-induced and thermally- generated degrees of freedom. 
This is illustrated with two qualitative points. First, at 
low fields within the ^-ordered state (Fig. 1), we can con- 
sider some effective "line-liquid" description (35). Right 
above T m (H), in the liquid phase, the self-consistent so- 
lution gives ($) = ($) l to be inserted in the effective 
interaction (p6[). Similarly, just below T m (H), we have 
($) = ($)s- In general, however, ^ (<&)s! This re- 
sult follows immediately from the gauge theory scenario 
(|]) since positional correlations of s vortices strongly in- 
fluence (<&). While both and ($)s are finite, as 
we cross T m (H), there is a discontinuous change in the 
average density and size of vortex loop and "overhang" 
excitations, resulting in different effective interaction ( J3q ) 
on two sides of the melting line (Fig. 4). The entropy 
jump at melting, AS", will receive significant contribution 
from such excitations. In fact, while AS — > as H — > 0, 
the critical fluctuations greatly enhance AS over the con- 
figurational entropy ,-pf field-induced lines. This provides 
natural explanation^] for the excess e ntrop y at T m (H) 
observed in low field thermodynamics.liirlia It should be 
stressed that this effect is different from the "microscopic 
entropy" contribution discussed by Hu and MacDonaldHl 
(see also Ref.Q). Such entropy arises from the electronic 
degrees of freedom and is reflected in the T-dependence 
of our GL coefficients (|l|) . This contribution is important 
in the high-field, LL regime. l3t3 In the 3D XY critical re- 
gion, however, such T-dependence is a minor effect since 
it involves T c and not the true transition temperature T c0 . 
The entropy contribution discussed here, arising from T- 
dependence of |(<E>)| 2 in Eq. (l36j) and ($) L ^ ($) s , is 



due to degrees of freedom of the superconducting order 
parameter itself and must be part of any proper theory 
of melting. 

Second, at higher fields, the self-consistent solution for 
Pg (or '!>) and <& leads to a rapid suppression of T$(H) 
far below T c q, and its subsequent merging with T m (H) 
for H > Hz (Fig. 1). This exposes a large region of 
the phase diagram where the Abrikosov vortex lattice 
melts directly into the true normal state. In this case 
= 0, while ($)s might still be close to its mean- 
field value. Such dramatic difference in the nature of 
these two phases, with the solid being rather unremark- 
ably mean-field like, and the liquid right above the melt- 
ing line exhibiting very strong fluctuations even at rather 
low fields, with vortex lines winding both along the field 
and in the perpendicular directions- evident in new 
simulations by Nguyen and SudbeO'EJ Note that such 
situation never arises in the mean-field based theories 
with only field-induced vortices,Q even after the appli- 
cation of our remedial procedure, since we would still 
have |($)s| 2 = |(<I>)l| 2 = m H =o\ 2 in Eq. ©. As 
argued in Sec. V, it does not appear possible to write 
an effective description of the true normal state in terms 
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of field-induced degrees of freedom only. For example, 
we could start again with our remedial procedure and 
argue that, even though it fails for q <C l/£, it still de- 
scribes the effective interaction of field-induced vortices 
for q ~ l/£, which is what matters most at T m (H). How- 
ever, within such a "line liquid" description, I do not see 
any simple way by which one could account for the part 
of AS associated with a discontinuous change in wind- 
ings across T m (H). In this region of higher fields, where 
T$(H) = T m (H), our gauge theory is becoming less and 
less "type-II" and amplitude fluctuations are becoming 
stronger in the GL theory, ft is likely that the low-field 
description discussed here and the LL theories now have 
an equal chance of providing reliable theory of the melt- 
ing transition: both routes, however, are certain to be 
most challenging. 
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APPENDIX A: DERIVATION OF THE GAUGE THEORY 

I present here a derivation of the gauge theory (Jfy. An abbreviated version is found in Ref.i. For simplicity, I 
consider the isotropic case, 7|| p,7_L = 7- 

Within the "helium model" ,c3 the partition function of the superconductor ([!]) can be written as: 



E II i II / 2*>ijexp(-f ); F v =7<l*| 2 > / d 3 r|V^(r)-^A| 2 + ]T £ / d Sl ^({ Xl Js L 

ATM w ' L=l ' w i w =l 



2 

where 



^EEEE /^/^^(ixd^-x^^^ , (37) 



VxA = H ; V-V^(r) = ; V x V<p(r) = 2tt V V / dxij(r - x, w [s i( J) . (38) 
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The partition function ( p7\ ) is a 3D counterpart of the familiar representation of the (continuum) 2D XY model in 
terms of its point-like topological excitations, vortices and anti- vortices. In 3D, the relevant excitations are loops/lines 
of vortices, classified by their global topology. Only vortex paths of unit vortex "c harg e" are considered since they 
are the important excitations in the critical region. The lattice regularization of (|37|) is a gas of non-intersecting 
oriented paths on a lattice which are either closed or can start/end only on sample surfaces. The lattice spacing is 
set by the characteristic "bending length" of vortex lines. Each individual step along a path takes a given energy 
to create and can be either up or down along x, y or z axis. These steps represent vortex segments and have a long 
range "directional" Coulomb interaction, operating only between the steps going along the same axis. This is the 
lattice version of the "Biot-Savart" interaction between vortex loops/lines of the continuum model. The background 
free energy, composed of the uniform condensation energy and "spin-wave" contributions, is not included explicitly. 

The summation in ( |37| ) runs over all distinct configurations of vortex line excitations of arbitrary length and 
shape. The index uj denotes different classes of oriented loops/lines which are distinguished by their global topology. 
For example, for periodic boundary conditions, when only closed loops are present, ui — (m x ; ; m^) , where 
m x,y z = 0, ±1,±2,... denote the winding numbers of a given loop around x,y,z directions. Finite closed loops 
correspond to m x = = mf = 0. Similarly, for free (periodic) boundary conditions along z (x, y) direction, 
u) = (m x ;m y ; 0, 0) denotes loops that wind in the x (y) direction while id = (0, 0; 0, 0; 0, 0) again denotes finite closed 
loops. In addition, there are vortex (antivortex) paths that traverse the system from z = to z = L z and "half-loops" 
which originate and terminate at the same z — or z = L z surface. JPx [u [sj u ] represents summation over all 
configurations of a given loop/line l u consistent with its global topology. Cjxi^fs;^]} in the first term of F v signifies 
that the integral of the gradient energy over the system excludes well-defined core regions associated with a given 
configuration of loops/lines. The second and third term represent core contributions: E^ = E c + ^({x;^ [s;J}) is a 

(2) 

"single-particle" term, with E c and Eb corresponding to the core line and bending energies, respectively. Vq { (|r — r'|) 
denotes "two-particle" effects of core overlap. These "two-particle" terms describe both the energy cost of core overlap 
and the entropic effects of keeping vortex and anti-vortex segments from annihilating each other. "Multi-particle" 
terms, arising from simultaneous overlap of more than two cores, can be neglected in the extreme type-II regime, where 
the average core size a is small compared to the average separation between vortex segments. Core contributions, like 
E c , Eb and Vg^(|r — r'|), can be computed in a specific microscopic model of vortex linesStii Their precise form 
is not needed for our present purposes since we will return to the GL representation at the end of this Appendix; it 

(2) 

suffices to know that E c and Eb are finite and the "interaction" Vq 1 (|r — r |) is short-ranged, of order a, and repulsive 

(2) 

on average. Without loss of generality, we could set Vq l / (|r — r'|) —> Vb,i<5(r — r ). 

With the uniform magnetic field present, the overall vortex "charge" neutrality demands that every configuration 
contains 7V$ field-induced vortex paths going from z = to z = L z . This fact is used to observe that the low 
temperature expansion of Z' (||) in terms of topological defects of the new order parameter <I>(r) and s vortices, is the 
same as ( |37| ) except for different prefactors: 1/(N$ + N a )\ in Z versus l/N$\N a l in Z' , where N a is the number of 
thermally-generated infinite vortex-antivortex paths which extend from z = to z — L z . This leads to a difference 
in entropy between two representations AS ~ Tln[(iV$ + N a )\/N$\N a \\. AS scales as L x L y , in contrast to the full 
entropy which goes as L x L y L z . Consequently, in the thermodynamic limit, AS does not affect the free energy per 
unit volume. In particular, within the ^-ordered state, Z (|l|) and Z' (^) have identical expansions: 

z = E -m*njn\ II f ztaok,]]! / £W Si ]exp(-^) ; 



j d 3 r\ V0(r) + U - ^A| 2 + ]T E c f ds h + ^ E c f ds 

Jc ' i =l i=l 



d 3 r / d 3 r'{[d (r)+d s (r)][ck(r')+d s (r')}V (2 \\r-r'\) + ^ . (39) 



Here {x; [s/ ]} denote finite closed loops of arbitrary length and shape. Half- loops attached to z = and z = L z 
surfaces are not included since they do not matter for the bulk thermodynamics. The symbol T>b indicates that the 
"single-particle" bending energy (^({x^ [s/J})) has been absorbed into the measure of the path integral. "Densities" 
and "currents" do, s and no jS are defined as: 

{d ,n }(r) = J2 [ dsi {l,^}5(r-xi [s lo ]) ; R,n s }(r) = ^ j d Si {l, ^-}5(r - Vi[ Si ]) , (40) 

Iq — 1 i— 1 
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and 



V-V0(r)=O , V x V<?i(r) = 27rn (r) ; V • U(r) = , V x U(r) = 27rn s (r) . (41) 

To understand the problem at low fields, we first consider the H = situation. The superconducting (Meissner) 
state is described by the H — * 0, iV$ — * limit of (p9|): 

oo JV (0) N m 

Z(H = 0)= J2 Jo)lII f^x io K](-^) ; F„(0) = 7<l*| 2 ) />r|V0(r)| 2 + J2 E c [ds lB + 

AT(0) =0 ' i = l C ?o=l 

+ i />r /dV^oWdoCr'^Gr-r'O + noW-noCrO^Cr-r'l)} . (42) 
The gradient term in F„(0) can be decoupled by a dual gauge field Ad(r), in the Coulomb gauge V • A c ; = 0: 

exp [_lM) jT d 3 r | W(r) |2] ^|p A(i (r)exp[|rf 3 r-(-mo-A d -^(Vx A d ) 2 )] , (43) 

where e 2 , = 87r 2 ^(|\E f | 2 ) /T is a dual charge. Eqs. ( [43| ) and ( fl2| ) have an appearance of a path integral over trajectories 
of relativistic charged quantum particles in a 2 + 1 dimensional Euclidean space [z being the imaginary time) coupled 
to the "electromagnetic" gauge potential A^ and interacting via short-range "density-density" and "current-current" 
interactions. Z(H = 0) describes the vacuum structure of such "electrodynamics" , with our vortex loops corresponding 
to worldlines of relativistic quantum bosons and describing virtual particle- antiparticle creation and annihilation 
processes in the vacuum. This similarity can be exploited further by using the particle-field duality to define the 
field-theory version!^ of Z(H = 0) (g§[4|): 

Jvy d (r) J VA d (r)cxp{- J d\[ml\V d \ 2 + \(V - tA d )V d \ 2 + ±g a \y d \ 4 + ^(V x A d ) 2 + ^-A 2 ]} . (44) 

^(r) is a field operator of these relativistic bosons. The gradient term in the action has been rescaled into dimension- 
less form so that short-range repulsion and dual charge assume their canonical dimensions: Vq5(y — r') — > <?o<5( r — r '), 
[go] — (length)" 1 ; e 2 d — > q d , [q 2 ] = (length) -1 . V\ has been dropped since it is irrelevant in the long wavelength 
limit. Its effect on critical behavior can be incorporated into the bare values of <?oi Q% an( i TO |>- Finally, the "bare 
mass" of Ad, M d , is absent in our problem (M d — 0). Finite M d reflects the presence of a gauge field minimally 
coupled to the original, superconducting order parameter. For example, if the condition k —> oo is relaxed and the 
real electromagnetic screening is restored in (|l|), M 2 — * fi e 2 /n, where e and ar e the real charge and magnetic 
permeability respectively. . .. . 

The expression (^) forms the basis for the "dual" picture of the 3D XY critical behavior.c§c2l In this picture, we 
are viewing vortices as primary objects and their field-operator, ^(r), as our order parameter, instead of the original 
^(r). We can think of *(r) in Ec. (0) as being the field-operator describing creation and annihilation of Cooper pairs. 
In the GL theory with H = (0), Cooper pairs have only short range interactions and it suffices to keep only the 
quartic term, describing the point-like repulsion, since the rest is irrelevant for the critical behavior. In contrast, the 
vortex excitations of \I/ interact via long range London-Biot-Savart forces, mediated by massless (M d = 0) A^ (|44|). 
Next, we can convert our neutral GL theory (Q) into the one with finite real charge, e, by introducing the fluctuating 
vector potential, A, as well as the electromagnetic field energy, (l/87r/zoe 2 )(V x A) 2 . Now, it is the Cooper pairs that 
have long range interactions, mediated by A, but the vortices in \& interact only through short range forces, due to 
the electromagneiic screening inducing a finite M d — fioe 2 /it in (fl4|). Precisely in 3D, there is an exact duality, at 
least on a latticejl3 between the form of this interaction for vortices and for Cooper pairsO This is what "inverted" 
stands for in the "inverted 3D XY" model. In the dual language of ^> d , it is the low temperature Meissner phase of 
the original superconductor that is symmetric ((^ d ) = 0), while the high temperature normal metal is the "broken 
symmetry" state of the dual theory ((^d) ^ 0). So, there is an inversion of the temperature axis. However, it is 
still the same symmetry, U(l), that is being broken and this implies that the thermodynamic exponent should be 
the same, v = v xy , both for e = and for e ^ 0. It is important to stress that the "inverted 3D XY" behavior of a 
charged-superfluid remains different from the 3D XY critical behavior of a neutral-superfluid, described by (Q) with 
H = 0, since they are associated with two different critical points. For example, the anomalous dimension exponent, 
of (^'(0)5 , *(r)), will be different in two cases, e = and e^0. 

Eq. (Q) describes the "true vacuum" state |0) of a Euclidean relativistic field theory. Particle (antiparticle) 
excitations are massive, m^, ~ 1/^ ~ 1/A*, where £ d is the dual correlation length associated with ^ d and 
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measures the average loop size. The actual value of m$ reflects both the cost in energy and gain in entropy arising 
from large thermally-generated vortex loops in the original problem. As we approach T c q from below, — > 0, and we 
enter the "false vacuum" state | /) of the theory ( f44| ) . Particle (antiparticle) excitations are now massless and infinite 
vortex paths proliferate across the system, as depicted in Fig. 2. This "false vacuum" is just the normal metallic 
state. If we introduce the particle (antiparticle) number operators, Np t A, |0) is an eigenstate of Np.A, Np\0) = 0|0) 
and Na\0) = 1 0} . In contrast, |/) is not an eigenstate of Np.A, and contains a finite average number of (anti)particles, 
(f\Np : A\f) — Np^a 7^ 0. Both |0) and |/), however, are eigenstates of the total vorticity operator Np — Na with 
eigenvalue 0, which insures Np — Na- For H small but finite, the ground state of (|4|) must still be an eigenstate 
of Np — Na but now the eigenvalue is N$. Starting from |0), such ground state |$) is naturally constructed by 
introducing N$ massive particles into the true, stable vacuum. We then have iVp|$) = A$|$) and Na\$) = 01$). On 
the other hand, starting from the "false vacuum" |/), the ground state \n) is formed by having additional N& particles 
added to an already present finite average number of massless particle- antiparticle pairs. \n) is not an eigenstate of 
Np t A and satisfies (n\Np t A\n) = Np.A ^ 0, where now Np — Na — Nq>. At T$(i7), a phase transition takes place 
between these two different types of ground state, |$) and \n), driven by the change in U(l) symmetry of the vortex 
system. 

We now return to Eq. (|3^). A finite density of s vortices produces two main effects on the loop "expansion" as we 
approach T$(iJ) from below. First, there is a long-range interaction between loops and s vortices which will influence 
the long-range correlations among the loops themselves. Second, there is a short-range effect of s vortices suppressing 
certain configurations of loops, through mutual contact interactions (an "excluded volume" effect). Intuitively, one 
expects the long-range effect to be essential for the critical behavior at low fields. Furthermore, the short range effect 
should be weak since the total number of vortex segments connected to s vortices forms a minority of all vortex 
segments. Based on these observations, we devise the following strategy: according to our basic assumption ii) and 
results of Sec. V, the system of s vortices below T$(H) can be viewed as equivalent, in the long wavelengthJimit, 
to an effective system of non-relativistic 2D quantum bosons in its superfluid state, in the sense of NelsonJa The 
collective modes of such a system are its "density" and "current" fluctuations, which is precisely how the loops couple 
to s vortices. The long wavelength effective action of these collective modes is described by two coupling constants, 
m s /n s and c 2 , where, in this boson analogy, n s , m s and c s are the superfluid density, mass and speed of "sound", 
respectively. We could compute these quantities explicitely, by starting sufficiently below T&(H). Here, however, we 
treat them as general parameters which characterize the long wavelength fluctuations of s vortices and whose ultimate 
values can be determined through their direct connection with the components of the helicity modulus tensor (Sec. 
V and Appendix B). We assume that the effect of mutual contact interaction between loop and s vortex subsystems, 
measured by some interaction strength V p , can be fully included into these ultimate values of m s /n a and c 2 , as well as 
into a renormalized loop- loop contact interaction, Vq and the loop core line energy E c . This amounts to reexpressing 



the short-range interaction in (39) as: 



d 3 r{Vo4(r) + 2V p 6d (r)6d s (r)+V pp d 2 s (r)} ; E c = E c + V Q (d s (r)) , E s = E c + V (d (r)) , (45) 



while the core line energies of loops and s vortices are also renormalized to include the average "excluded volume" 
effect. E s and V pp will shortly be subsumed into m s /n s and c 2 . We set V p — > in ([l5]) and proceed to derive the long 
distance description of (|39|). Indeed, we find that the long-range interactions lead to a major change in the behavior 
of the system once H is finite: most importantly, the long-range interactions between vortex loops are "screened" by 
fluctuations of s vortex lines. The amount of "screening" is determined by m s /n s and c 2 . Once the critical behavior 
associated with this "screening" has been understood, we reinstate the residual contact interaction between loop and 
s vortex subsystems and test for consistency. We find, both within the e expansion and perturbative RG in fixed 
dimension D — 3, that such residual coupling is irrelevant for 3 < D < 4, i.e., it does not lead to any new relevant 
terms in the effective action, apart from those alread y p resent. Our procedure is therefore self-consistent. 

The first step is to decouple the gradient term in (|3S|) by using Ad, except that now no — > no + An s in Eq. (|43| ) 
(An s is defined above Eq. (IJ)). The s vortex part of (J39|) becomes: 

/ ZW*]exp{-5^ / ds 4 -^ [ d\d 2 s {v)~ l f d 3 rAn s (r)-A d (r)} . (46) 

This part is now reexpressed in terms of s vortex "density" and "current" fluctuations. Sufficiently below T$(iJl the 
overhangs are small and s vortex lines are almost directed. We can use a non-relativistic boson analogy of NelsonS and 
replace J dsi m f^* dz(l + ^(^) 2 ), where n(z) = (xi(z), yi(z)), while absorbing the effect of small overhangs into 
the effective mass, m s ~ E s /T. Similarly, d s (r) — > n s ii(r) — * ^(z)) and n s j_(r) — * Y^ii^il dz)S(f '— fi(z)). 
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We then introduce a field operator of these non-relativistic bosons, 4' s (r, z), replace the path integral in flf6| ) with 
the functional integral over ^ s (periodic— <pr free boundary conditions give the same result in the thermodynamic 
limit), and follow the standard procedurecS for deriving the "hydrodynamic" action of s vorticesS the phase <p a (f, z) 
and amplitude 7r s (r, z) are introduced via ^> s — \*f? s \ exp(iip s ) and tt s = \^ s \ 2 — (\^s\ 2 ), where (|^s| 2 ) is evaluated 
self-consistently. After integration over w s , the "hydrodynamic" action for the phase becomes: 

^M^(V \\<Ps -A d|| ) 2 + -^(V^ S -A d± ) 2 + (-..)} , (47) 

where (• ■ •) denotes higher powers of Vy s — A d and higher order derivatives. The "mean-field" part of the action is 
absorbed into the background. The action ( |47| ) is decoupled using real fields 8p{v) and j( r ) = (jx,jy)' 

Jv5p(r)Jvj(r)cxp{Jd 3 r[-i8pA dll -ij-A d ^+i8pV^ s +ij-V^ s ~ 1 ^5p 2 -^j 2 ]} , (48) 

which finally results in Eq. ( p9f ) expressed in terms of the true collective modes of the superfluid s vortex system, its 
"density" 8p and "current" j fluctuations: 

oo N (0) N m ~ 

Z ^ E i II fVwM I ™p(*) [ "fa) [ VA d (r)eMS) ; S= £ [ ds lo + 

AT(0) =0 ' l = l J J J J Jo=1 J 

+ J d\{^d 2 (v) + m • A d + i5 P A d || + ij- A d± + ^5p 2 + ^-j 2 + W(Sp,j) + l(Vx A rf ) 2 } , (49) 

where the functional integral must be appended by the continuity condition, V||<5p + Vj_ ■ j = 0, which follows from the 

integration over tp s in (ff8|). W arises from (• • •) terms in and contains powers higher than quadratic in Sp and j, 
as well as assorted derivatives. In particular, odd powers of Sp are present, like ~ Sp 3 , reflecting the An g ii — ► — An g ii 
asymmetry of the original problem (|39|). Such higher order asymmetric terms are essential for description of the 
melting transition but, as shown below, are irrelevant at the ^-transition, provided the latter is continuous. Note that 
the above relatively simple dependence of ( [49|) on Sp and j holds only at distances > t. This is precisely what we are 
interested in as we approach T$ (H) . Eq. fl49| ) captures an essential effect of a finite field on the loop "expansion" : 
fluctuations of field-induced s vortex lines result in the "screening" of the "Biot-Savart" interaction between the loops. 
This "screening" is manifested by A^ gaining a finite "mass", M 2 ~ n s /m s c 2 , M\ ~ n s /m s , after integration over 

Sp and j. The effect of finite magnetic field in the original problem (g) is now stored in the finite values of M| and 
Mj^. As one attempts to create a small density of very large loops in (f49j), upon approaching T$(H) from below, their 
effective line tension and mutual interactions will be essentially influenced by such "screening" . 

The importance of this "screening" mechanism is particularly apparent in the dual representation (^) . The action 
(E^) becomes: 

f d\[ m % |$ d | 2 + |(V - iA d )$ d \ 2 + ||$ d | 4 + i6 P A dz + ill d + ^Sp 2 + ^-f + W|i(Vx A d ) 2 ] }, (50) 

where the meaning of $d(r), m$, and go is evident in light of the discussion surrounding Eq. (^4|). For simplicity, I 
suppress the anisotropy in the second (gradient) term of (^o|) which generically arises (even if 711 = 7_i_ in Q) from 

the interaction of loops with s vortices. The finite mass of A^, generated by the integration over Sp and j, "cuts off' 
the long-range "Biot-Savart" interactions present in zero field (|44|). This "screening" causes decoupling of the dual 
gauge field and transformes the critical behawior from a "charged" (|44| ) to a "neutral" (|m]) dual superfluid. This is 
the inverted (anisotropic) 3D XY behaviorHca of the dual theory, hinting at the presence of a massless gauge field in 
the original "superconducting" formulation (Q), as discussed below Eq. ( ff4|). 

To investigate the critical behavior of ( |49"| ) in more detail, we generalize (|50|) to arbitrary dimension D. Simultane- 
ously, we restore in the action the residual part V p (|4^) of the contact interaction between loops and s vortex lines, 
i.e., the part not already incorporated into the values of coupling constants appearing in (p0|): 

J d D r[m% \<S> d \ 2 + \{d„ - iA p )<S> d \ 2 + igol^l 4 + g p \<S> d \ 2 J + U^A" + ^^J" + + ^2 Vl . (51) 

Here g p denotes V p rescaled to its canonical dimension, p = 0, 1, 2, . . . , D — 1, W( J p ) is the generalization of W(Sp,j), 
and F pv — d p A v ~ d v A p . The functional integration runs over fields & d , A p and J p and includes a constraint 
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d^J^ = 0. The integration over J p generates a finite mass M M= o == Mt\, M M ^o = M± for the dual gauge field A p . 
The theory, however, retains gauge invariance implying that the combination qdA p must be an RG invariant. This 
in turn sets the canonical dimension of g p to [g p ] — (length) D ~ 3 . Therefore, the term g p \$d\ 2 Jo is irrelevant within 
e expansion around the upper critical dimension D — 4. Similarly, all higher powers and derivatives of J p appearing 
in W, are irrelevant as well. For example, the canonical dimension of the Jq coupling constant is (length) 2 - 0-3 . The 
relevant couplings below D = 4 are go, q\ and Mim. Due to finite M»±, A p decouples and the j3 function for go is 
the same as that of the neutral (gy = 0) complex $ 4 -theory. In this case the e expansion is expected to hold down 
to, and include, D — 3, where the critical behavior of the <I> 4 -theory should be that of a (inverted) 3D XY model, in 
agreement with our earlier assertion. The same conclusions can also be reached within the perturbative RG in fixed 
dimension D = 3. Here g p is marginal at the "engineering" level, the gauge field again decouples due to finite Af|i j_, 
and we can compute the relevant (3 functions at the one-loop order and to the leading order in g p : 

Pa(9o, 9 P ) = dl Qg (p) = So + Ciffo ; Pp(9o, g P ) = dl pg(p) = C 29og P , (52) 

where go,p(p) are the dimensionless running coupling constants and are (regularization dependent) numerical 
constants which are both positive, C12 > 0. At the (inverted) 3D XY critical point go = 1/Ci and therefore j3 p > 0, 
indicating stability of our assumed g p = fixed point against residual g p perturbation. The above results allow us 
to conclude that the critical theory (EG) (with W = 0) remains valid and that the effects of V p can be included by a 
proper choice of relevant couplings, as originally assumed. The presence of long-range interactions between vortices, 
mediated by A<j, is essential for the validity of this argument. 

One step remains: as T — > T$(H), some overhangs attached to s vortex lines become very large and we might 
doubt the accuracy of the straightforward non-relativistic boson analogy approximations below Eq. (^). However, 
throughout the ^-ordered state, the s vortices remain "massive" and there should always exist a suitably defined 
quantum system of non-relativistic 2D bosons whose long "distance" (x, y) and "imaginary time" (z) behavior faith- 
fully emulates that of s vortices. We therefore expect that the overall symmetry of ( (49[ ) remains preserved at T$>(H). 
This leads to the generalization of (|^j): 

oo N (0) N m ~ 

Z W - E ill <fv bXl0 [s l0 ] f VV(r) f 2?A d (r)exp(-S) ; S = ]T -f f ds lo + 

rfM^o(r)+m .A d + i V.A d + ^Vf + ^ivi + ^(VxA (i ) 2 } , (53) 

where V(r) describes long distance (3> tj fluctuations of s vortex "currents" n s (r) ( p0| , p] ) and satisfies V • V = 0. 
K\\j_/T now play the role of m s c 2 s /n s and m s /n s in ( |4S| ) and fully include the effect of overhang configurations as 

T — > Tg,(H). At present, we cannot compute Kn ±(T, H) (nor E c and Vq) from first principles. This would require 
an analytic solution to the problem of large overhangs, something far beyond the scope of this paper. However, 
if we start with the general form (p3|), we can determine various parameters appearing there by connecting them 
self-consistently to directly (numerically or experimentally) measurable physical quantities. For example, K n ±(T, H) 
can be extracted from the components of the helicity modulus tensor (Appendix B) or the fluctuation conductivity 
(Sec. VI). We should therefore consider ( [531 ) a self-consistent, perturbative RG description of the ^-transition. 

We can now enforce the constraint V • V = by introducing a gauge field S(r): 2irV -> V x S. V ■ S = 0. 
Alternatively, we can integrate over V, obtain the mass term for A^, and then decouple it by introducing S. The 
final result is: 

oo N i0) N m ~ 

Z{ ~ H ) - E ill <f>V b xM [ ^S(r) f 2?A d (r)exp(-S) ; S = ]T [ ds lo + 

dM^dfc) + in ■ A d + ±(y x S) • A d + §(V x S)f + |^(V xS)l + i(Vx A d ) 2 } . (54) 

This is just the vortex loop expansion ([l||,[l3|) of the Meissner phase of a "superconductor" described by an order 
parameter $(r) and coupled to the gauge field S (V • S = 0). The GL functional of such a superconductor is 

F cS = 5|<i>| 2 + %\{V P + *S M )<i>| 2 + fl*! 4 + ^(V xS)j + |(Vx S)i , (55) 



24 



which is precisely the gauge theory (Q). a, $ and 7 M are some suitably renormalized GL coefficients which can be 
determined phenomenologically. Note that we have now restored the anisotropy in 7^, arising both from the bare 
anisotropy (7|| 7^ 71 in Eq. (|])) and the anisotropy induced by the interaction of loops with the s vortex background. 



APPENDIX B: GAUGE THEORY AND THE HELICITY MODULUS 

Here we consider the connection between K± t \\ appearing in Eq. M)_and the helicity modulus tensor T(q). The 
conventional definition of the components of T(q) can be found in Ref.E3: 



T„„(q) = V 



5 2 F 



<5a„(q)5a M (-q) 



(56) 



All quantities appearing in ([56]) are defined in Sec. V, below Eq. (|T^). We will limit our attention to the isotropic 
case (71 = 7|| in Eq. (|l|)) whenever we consider the H ^ situation. The anisotropic case with finite field requires far 
more extensive algebra and can be reconstructed by combining the discussion below with illuminating presentation 
in RefM 

We first evaluate T(q) from the original GL theory (P and start with the H = (A = 0) case. We consider the 
situation right above T c o, so there is no superconducting long range order. After adding small a to the gradient term, 
we can expand the free energy up to second order in a: 



F[V x a] = F[0] 



d 3 r 



2c? 



X± (V x a) 



2?Xy(Vxa)| 



(57) 



This equation requires a brief explanation: e is the real electric charge and c is the real speed of light, appearing in 
(M. F[0] is just the original free energy of the GL theory with H = and a = 0. The full free energy with small 
a(r) is written as a functional of V x a only. This is required by the gauge invariance of ([l]). Two additional terms, 
proportional to the perpendicular and parallel components of (V x a) 2 , represent the leading corrections in powers 
and derivatives of V x a. The subleading contributions are denoted by dots. These subleading terms are unimportant 
in the long wavelength limit. 

Xj_ and x\\ are some functions of T and are different from each other in the anisotropic case, 7|| j±, while 
Xi_ = Xll = x if the superconductor is isotropic. As T approaches T c q from above, these functions take the following 
form: 



CT6 



CT 



1 



(58) 



= Co_L,||^ v are the superconducting correlation lengths and C is an unknown universal constant, intrinsic to the 
GL theory ([l]). Note that e 2 x±,\\/c 2 is just the perpendicular (parallel) magnetic susceptibility. 

The components of the helicity modulus tensor in the long wavelength limit (q — > 0) are uniquely determined by 
X_l,||. Conversely, the measurement of the long wavelength "tilt" and "compression" helicity modul£9 determines 
X_l ||- In general, from definition (pq), we find: 



rY M „(q) = X^pa^pPuQaQp 



(59) 



for the isotropic case, while for the anisotropic situation x_l X\\ : 



— T /^(q) = (X|| ~ X^)^za^z(3yq a qp + X^pa^pf3^q a q/3 ■ 



(60) 



e Q /3 7 is the Levi-Civita symbol and summation over repeated indices is understood. e 2 /c 2 is factored out for later 
convenience. 

After this preliminary discussion, we go to the case of interest, finite H in Eq. ([l]), and limit our consideration to 
the isotropic situation 7^ = 7|| . We introduce small a into Eq. ([!]) and expand to second order in V x a: 



F[H,V x a] - F[H,0] + \ d 3 r ^X±(V x a) 2 



e 

2c 2 



X||(Vxa) 2 



+ 



(61) 



25 



where F[H, 0] is now the free energy of the GL theory ([j]) at finite field H and a = 0. We are again focusing on the 
"normal" , i.e. not superconducting state, in accordance with assumption ii) of Sec. II. The above expression looks 
very much like j57| ) but there are following significant differences: the expansion is anisotropic, x_l 7^ Xh even though 
our superconductor is isotropic. The reason for this is finite field H along the z-axis which reduces spherical symmetry 
of the H = situation down to cylindrical. The finite field also breaks the "up-down" symmetry along the z-axis. 
This is manifested by the subleading corrections, denoted by dots in d6l]), containing in general odd powers of V x a 
(the leading such term is cubic); such terms were prohibited by symmetry in the H = case. Again, by combining 
definition (|56| ) and Eq. (|6l|), we arrive at the expression for the long wavelength limit of the helicity modulus: 

c 2 

(q) = (x|| - X±.) £ z a ^z/3uq a q0 + X±e pa ^pf3»<laq0 ■ (62) 

The components of the helicity modulus tensor for the finite field (isotropic) case are determined by X±,\\ which are 
some functions of T and H . It is tempting to conclude that: 

X ±^CT^(T,H) , XW^CT ^pJ^ , (63) 

where £_u|(T, H) are now superconducting correlation lengths at finite field. This result is plausible on physical 
grounds, expressing the fact that, with H =/= 0, the superconducting correlation lengths are now finite in the "liquid" 
phase, limited by magnetic length I and consequently the helicity moduli vanish in the q — > limit. We are simply 
making the assumption that the same length that limits the range of superconducting correlations appears in the 
coefficient of the g 2 -term in the helicity modulus; this assumption is known to be correct for the H = case (|5g|). 
Unfortunately, I am unable to provide a mathematical proof that the conjectured result ([33j) is exact. Instead, I define 
perpendicular and parallel "screening" lengths A±(T,H) and Ay(T, H) by: 

X±=CTA ]] (T,H) , X|| =CT ^(r'g) • (64) 

Note that there is a one-to-one correspondence between these "screening" lengths Ami and x_i_,|| and, in turn, between 
|| and the long wavelength helicity moduli ([62]). Since Aj_j| are purely thermodynamic quantities they satisfy scaling 
laws of Sec. IV, just like the superconducting correlation lengths: 

A_ L ,||(T,JT)=0^ i ||(gg) , tL,\\{T,H)=m 6 ±ii (4) , (65) 

where "dimensionless charge" q% = f (T, H = 0)/£ oc E x l % j\t\ v is our scaling variable of Eq. @ and ^ and TZ^ ,, 
are the screening length and correlation length scaling functions, respectively. In the H — > (q$ — > 0) limit, all these 
scaling functions, TZ± ^ (<7q) and 1Z^_ ,, (qq), go as g 2 . We now make the following assumption: around T$(H), the ratios 

A±(T,g) Kj{gl) AjigVg) ^(go 2 ) 

are some unremarkable smooth functions of the scaling variable g 2 . In particular, the non-analytic drop in the 
coefficient of the g 2 -term in the helicity modulus, which takes place as we cross the T^{H) transition line (|2l]) and 
which is directly reflected as a non-analytic decrease in the screening lengths A±\\(T,H), is manifested also in the 
superconducting correlation lengths £j_ 11 (T, H). This assumption, which appears justified physical grounds, was used 
in the discussion of fluctuation conductivity (Sec. VI). 

Finally, we are in position to discuss our anisotropic gauge theory of Eq. (Ej) . Small a added to the external vector 
potential A in the original GL theory (|l|) translates into a small vector potential (e/c)a added to our fictitious gauge 
fieldS in Eq. (|). Since we are integrating over S, it is useful to define new gauge field S n = S + (e/c)a and integrate 
over S n in the partition function. The effect of this is to move (e/c)a from the covariant gradient terms, |-D M $| 2 , to 
the "gauge field energy" A'^V x (S — (e/c)a)] 2 in Eq. (|j). We now expand the free energy of the gauge theory, F e g, 
to second order in V x a, following the same philosophy as in Eq. (|6f|): 

r 2 2 

F eff [e x ,e||,Vxa] =F eff [e ±) g||,0]+ / d 3 r[^/C ± (V x a)i + ^/C||(V x a)jj] + . . . , (67) 

where 
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£_L,|| = - 



-ifr- lim / d 3 (r-r') e iq ' (r - r,) ((V x S(r)) ±i y(V x S(r')) ± ,||) . (68) 



The anisotropic charges e± i \\(T,H) and coupling constants K± I \\(T,H) are defined in Eqs. and (|4|), respectively. 
The thermal average (• • •) is over the gauge theory defined by the free energy functional !F e ff (^). Combining (|56[ ) and 
(|68|), we get, as before: 

c 2 

^^^(q) = ~ &l-)tza^z[3vquqf3 + K, 1.Z pa.^ P Pv<la% ■ (69) 



Comparing this to the general expression for the helicity modulus of the GL theory at finite field, given by Eq. (|6z 
we conclude that xx,|| — /C_|_ ) 11 - 

A very important point: below T$(H) we have K±_ 11 = id ||- This is a mathematical consequence of the following 
physical picture (Sec. V). In the ^-ordered state, only the field-induced (s) vortex lines can "screen" the test vector 
potential a(r) . All thermally generated vortex loops are of finite size and cannot contribute to the "screening" in 
the long wavelength limit, described by /Cx,|| m Eq. (|67|). This implies that the new order parameter $ is finite 
and therefore ((V x S)^ i.) ( |68|) vanishes in the long wavelength limit, as ~ q 2 /\(&)\ 2 . Therefore, we can uniquely fix 
the coupling constants K± \\(T,H) (and corresponding anisotropic charges e±u(T,H)) that enter the gauge theory 
description (||), by connecting them directly to the components of the helicity modulus tensor right below T$(H) (or, 
more precisely, for T$,(H) — T — > + ). In particular, the functions Cx,\\(T,H), introduced below Eq. (|J) and referred 
to at various points in the main text, follow from the general expression (|6~i|): 

c x (T,H)=CHf(ql) , C ||(T,g)=c M^ , (70) 

where 72.^; m(3o) are the scaling functions introduced in Eq. ( |65|) and are to be evaluated below Tg,(H). The gauge 
theory scenario predicts that the fundamental anisotropy ratio cii/cj_ takes on a universal value along T$(H). 

Above T$(i?), in the true normal state, infinite vortex loops proliferate across the system and can contribute to 
screening. ((V x S)^ ,,) in Eq. ( |68| ) becomes finite and ~ This causes the non-analytic drop in the g 2 -term of 



the helicity modulus and the corresponding screening lengths (J64|), just as discussed in Sec. V (see Eq. (21)). 

The reader should note that the set of results presented in this Appendix, connecting the properties of the gauge 
theory description (Q) around T$ (H) to the general long distance form of the helicity modulus and screening lengths 
of the GL theory at finite field (|l|), is not only physically transparent and appealing but also exact, provided our 
assumptions i)-iv) (Sec. II), are satisfied. 
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FIG. 1. Proposed H — T phase diagram for the critical region of extreme type-II superconductors. The dashed region 
denotes a crossover from the Gaussian regime, where amplitude fluctuations are strong, to the critical 3D XY-like regime, 
where amplitude fluctuations are suppressed. Within the 3D XY-like critical regime, the London-type vortex loops and lines 
with tight cores are well-defined excitations. Along the temperature axis, this critical region is bounded by the mean-field T c . 
Along the field axis, the critical region is boundadibv H s (~ Ht) (|^). Above H a , the physics of GL theory (Q) is dominated by 
the formation of Landau levels for Cooper pairs.L3'E3 The "^-transition, or the vortex-loop "expansion" transition, T<s>(H), and 
the vortex lattice melting, T m (H), occur simultaneously for H > Hz ~ H s . This is a first-order transition from the Abrikosov 
vortex lattice (VS) directly to the normal state (N). Below Hz, T<s>(H) and T m (H) split into two separate transitions and merge 
again only at the true zero-field superconducting transition, T c o, as H — > 0. For H <C H B , the transition at T<s,(H) is continuous, 
while the vortex lattice melting transition remains first-order. The intermediate phase ($), below T<s>(H) but above T m (H), is 
not a superconductor = 0), but it differs from the true normal state (N) by a new type of long range order, characterized 
by the "superconducting" order parameter $(r) (|^,^). Only N$ field-induced vortex lines traverse the system along the field 
direction in this $-ordered state, while the average size of thermally-generated vortex loops is finite. In the true normal state 
(N), the $-order is destroyed as numerous additional vortex paths "expand" across the system in all directions. 



FIG. 2. A schematic representation of the H = transition in an extreme type-II superconductor. The low temperature 
Meissner phase (T < T c o) contains only finite vortex loops. In the high temperature normal state (T > T c o) these loops connect 
and "expand" across the system, leading to loss of phase coherence and finite dissipation. Two forms of superconducting order, 
described by \I> and $, are equivalent here. For clarity, the vortex paths are drawn far smoother than they actually are near 



FIG. 3. Characteristic configurations of the system a) below and b) above T<s>(H) (but always above T m (H)). a) Field-induced 
vortices (depicted in blue) wind all the way across the system along the field direction but only undergo effective "diffusion" 
in the transverse direction, b) After the loop "expansion" at T<s,(H) this effective transverse "diffusion" is destroyed, as 
field-induced vortices can "hitch" a ride all the way across the system in the xy-plane by connecting to thermally-generated 
infinite loops present in the true normal state. Note the presence of "vortex tachyons" (depicted in red) which wind only in 
the xy-plane. Again, for illustrative purposes, the vortex paths are drawn far smoother than they actually are near T<s>(H). 



FIG. 4. Characteristic configurations of the system a) below and b) above T m (H) (but always below T<s>(H)). a) Field-induced 
vortices (blue) execute small oscillations around their equilibrium positions. Thermally-generated vortex-loops (red) are small 
and rare, b) Above the melting transition thermally-generated loops discontinuously grow larger and more numerous, although 
they still remain of finite size. This discontinuous change in the state of the loops accounts for (3>)s 7^ at the simplest 

mean-field level. 
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